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ABSTRACT 

A  number  of  simple  models  of  diffusion  electrodes  are  examined 
theoretically.  The  results  are  expressed  in  forms  whereby  it  should 
be  possible  to  test  the  models  experimentally,  and  make  predictions 
regarding  the  effects  of  variation  of  the  pertinent  parameters. 

It  is  not  expected  that  the  theories  will  completely  explain 
operating  data  on  complex  electrode  structures.  It  should  be  possible 
to  test  them,  however,  on  experimental  equipment  designed  to  approximate 
to  the  models  chosen. 


(1)  INTRODUCTION 


Although  primary  and  secondary  cells  have  employed  porous  electrodes  for 
many  years,  there  appears  to  have  been  little  work  on  describing  the  exact  Influence 
of  the  porosity .  Presumably  this  is  due  to  the  complex  nature  of  the  electrochemical 
reactions  in  such  cells  and  the  difficulty  of  describing  paste  electrodes.  The 
recent  surge  of  interest  in  fuel  cells  has  led  to  Increased  study  of  porous 
electrodes,  with  some  hopes  of  applying  reasonably  exact  mathematical  treatments. 
Fuel  cell  electrodes  are  usually  hard  and  solid  and  retain  their  configuration 
during  operation,  since  they  serve  to  provide  catalytic  reaction  surface  rather 
than  to  provide  storage  of  reactant.  Justi  £t  (1)  have  presented  a  limited 
analysis  of  porous  fuel  cell  electrodes.  A  similar  analysis  was  performed  inde¬ 
pendently  by  Austin  (2).  This  analysis  is  incorporated  into  this  report  since  it 
forms  the  basis  for  more  detailed  treatments.  Stender  and  Ksenzhak  (3)  have 
discussed  the  limiting  case  where  polarization  is  due  solely  to  ohmic  effects 
within  the  electrolyte, ift  the  pores  of  the  electrode  and  no  depletion  of  reactant 
occurs  in  the  electrolyte.  We  will  show  that  this  case  is  only  rarely  appli¬ 
cable.  The  transient  case,  which  is  difficult  to  analyze  when  appreciable  amoi-its 
of  reactant  are  adsorbed  on  the  surface,  has  been  crudely  treated  by  Uhrbach  (4) 

A  number  of  other  workers  (5-12)  have  presented  studies  of  various  porous  systems, 
but  since  the  work  is  as  yet  available  only  in  abstract  form,  it  is  not  possible 
to  critically  analyze  their  findings  (although  reference  to  their  presentations 
will  be  made  in  this  report) . 

One  of  the  difficulties  in  describing  the  mode  of  operation  of  diffusion 
electrodes  is  that  the  complexity  of  the  systems  leads  to  final  solutions  which 
contain  several  parameters  in  "awkward"  algebraic  forms .  Thus  it  is  conceivable 
that  by  suitable  selection  of  the  values  of  the  parameters,  several  different 
theories  can  be  made  to  fit  a  given  set  of  experimental  results.  In  order  to 
test  certain  assumptions,  it  is  valuable  to  consider  simple  redox  reactions  at 
porous  electrodes;  the  reactions  have  fairly  well  defined  kinetics  and  a  two  phase 
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system  results.  Theories  can  be  relatively  easily  checked  for  such  cases. 

:  The  concepts  can  then  be  extended  to  three  phase  gas  diffusion  electrodes 

with  more  certainty  after  they  have  been  proved  correct  for  the  simpler 
’  systems. 

In  general,  we  are  looking  for  a  model  (or  models)  of  a  porous 
electrode  which  will  satisfy  two  requirements.  It  must  be  sufficiently 
rigorous  to  explain  quantitatively  experimental  results  from  a  simple 
system,  for  example,  a  single  pore  of  known  radius  and  length.  In 
addition,  it  must  be  general  enough  to  explain  qualitatively  the  behaviour 
of  complex  systems  such  as  sintered,  double  layer  structure,  Raney  electrodes 
(13)  or  carbon  electrodes.  It  is,  of  course,  possible  to  write  differential 
equations  covering  every  conceivable  factor  affecting  the  results.  It  is 
also  possible,  in  theory,  to  solve  the  equatipns  numerically  u61ng  a  digital 
computer.  But  we  are  likely  to  end  up  with  solutions  which  have  so  many 
unknown  parameters  that  they  can  be  used  to  fit  any  set  of  experimental 
results.  It  is  my  opinion  that  the  better  way  to  approach  the  problem  is 
to  start  with  the  simplest  possible  models,  solve  the  equations,  and  perform 
experimental  work  specifically  aimed  at  testing  the  models.  By  examining 
the  problem  at  several  simplified  limiting  conditions,  it  is  often  possible 
to  show  that  a  particular  effect  is  usually  negligible.  It  is  also  my 
opinion  that  fairly  lengthy  explanations  of  the  methods  of  setting  up  the 
differential  equations  are  not  superfluous.  The  equations  and  their 
solutions  may  be  valuable  to  people  who  have  not  had  much  experience  in 
writing  down  these  equations,  and  the  extra  space  used  in  describing  the 
meaning  of  the  equations  is  not  wasted. 

(2)  PHYSICAL  DESCRIPTION  OF  ELECTRODES 

Several  simplifying  assumptions  must  be  made  about  the  structure 
of  porous  electrodes  before  it  is  possible  to  apply  simple  equations.  These 
assumptions  may  not  be  completely  true,  but  they  are  necessary  to  form  a 
model  amenable  to  treatment.  For  example,  it  is  not  generally  possible  to 
get  a  porous  electrode  which  is  perfectly  homogenous  in  pore  size.  A 
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small  pore  will  have  a  greater  surface  area  to  cross-sectional  area  ratio 
than  a  larger  pore,  If  both  pores  are  held  at  the  same  polarization  so 
that  current  per  unit  area  is  the  same,  and  reactant  allowed  to  diffuse 
in  at  one  end  of  the  pore,  the  smaller  pore  must  diffuse  in  a  proportional¬ 
ly  greater  amount  of  reactant  per  unit  cross-sectional  area  than  the  large 
pore.  Consequently,  the  concentration  of  reactant  will  fall  off  more 
rapidly  with  distance  from  the  entrance  mouth  for  the  small  pore.  A 
bundle  of  parallel  tubes  of  varying  diameter  could  be  used  as  a  porous 
electrode.  Theoretical  equations  for  current  density  versus  polarization 
for  each  pore  could  be  obtained  and  using  the  common  boundary  conditions 
of  uniform  concentration  and  polarization  at  the  pore  mouths,  the  sum 
effect  could  be  found,  providing  the  size/number  distribution  was  known. 

Ohmic  potential  gradients  and  concentration  gradients  \(within  the  electro¬ 
lyte  in  the  pores)  would  vary  with  each  pore.  However,  if  one  examines 
a  porous  carbon  or  sintered  metal  electrode  under  a  microscope  it  is 
apparent  that  this  model  is  not  adequate.  These  structures  are  spongelike, 
with  the  pore  spaces  interlinked  in  many  directions.  Any  variation  of 
concentration  gradient,  or  electrolyte  potential  gradient,  existing  along 
two  different  size  pores  will  tend  to  even  out  by  lateral  mass  transport 
through  the  Inter-llnklng  connections  of  the  pores.  This  discussion  leads 
to  the  first  assumption.  It  will  be  assumed  that  the  pores  are  so  highly 
inter-linked  that  the  potential  and  concentration  at  any  given  penetration  into 
the  electrode  are  constant.  This  is  probably  a  good  approximation  for  well 
constructed  electrodes.  An  electrode  with  a  crack  (or  large  pin-holes) 
through  it,  on  the  other  hand,  will  behave  somewhere  in  between  the  extremes 
of  completely  interlinked  pores  and  discrete  tube  bundles. 

Another  assumption  is  that,  in  general,  a  porous  electrode  has  a  macro - 
porosity  and  a  microporosity.  The  macroporosity  is  the  space  between  particles 
or  material  which  can  be  observed  under  an  optical  microscope.  Microporosity 
is  very  fine  porosity  existing  within  or  on  the  surface  of  the  large  material. 
For  example,  in  a  baked  porous  electrode  carbon,  the  microporosity  is  made  up 
of  Inter.fcry^telllte'  voids  in  the  filler  and  fine  channels  or  bubbles  in  the 
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carbonized  binder.  For  a  platinized  platinum  electrode,  the  voids  between 
the  base  platinum  particles  form  macroporosity,  while  the  fine  surface 
porosity  of  the  platinized  layer  is  the  microporosity.  In  Raney  metal 
electrodes,  the  microporosity  is  the  fine  pore  structure  within  the 
particles,  produced  by  dissolution  of  the  soluble  alloying  component. 

The  macroporosity  can  be  considered  as  having  a  definite  cross-sectional 
area  and  geometric  perimeter  at  any  cross-section  of  the  electrode.  Normal¬ 
ly,  the  microporosity  represents  only  a  few  per  cent  or  less  of  the  total 
empty  volume  and  we  may  equate  macroporosity  to  G,  the  total  porosity.  The 
total  perimeter  area  of  the  macroporosity  represents  the  area  through  which 
mass  transport  to  the  reacting  surfaces  can  occur,  and  may  be  denoted  by  A 
cm^  per  cra^  of  electrode.  The  total  area  for  reaction,  S  cm^  per  cm^  say, 
can  be  considered  as  A  times  a  roughness  factor  f, 

S  =  Af  cm^  cm  ^  (1) 

f  is  all  the  area  associated  with  unit  geometric  macroporous  area  A.  For 
completely  smooth  particles  (or  tubes),  with  no  microporosity,  f  =  1.  For 
a  particle  with  highly  developed  interior  or  surface  area,  f  will  be  large. 

Another  assumption  made  is  that  the  conductivity  of  the  electrode 
material  is  large  compared  to  the  conductivity  in  the  electrolyte.  For 
small  experimental  electrodes,  ohmic  effects  within  the  electrode  material 
will  be  neglected. 

It  will  also  be  assumed  that  the  electrodes  investigated  will  be  plane 
electrodes,  with  uniform  current  density  in  the  external  electrolyte  and 
negligible  edge  effects. 

(3)  REDOX  SYSTEMS  (IWO  PHASE  SYSTEMS) 

3.1.  General  Discussion 

In  the  consideration  of  redox  reactions  at  a  porous  electrode 
the  problem  is  simplified  since  we  need  consider  only  transport  in  the 
electrolyte  phase.  It  is  assumed  that  any  surface  intermediate  is  present 
in  small  quantities  (9  -»  0,  where  9  is  the  fractional  surface  coverage) . 

The  activity  of  bare  surface,  1-0,  is,  therefore,  close  to  1  and  can  be  taken 
as  constant  over  a  wide  voltage  range.  The  rate  equation  can  be  expressed  as 
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(14), 


_  On.FTl/RT  o  -(l-a)n,FTl/RT 

-17'  .  ) 

i  is  current  density,  i^  is  the  exchange  current  density  for  concentrations 
(activities)  of  reactant  and  product  of  R^,  P^,  respectively.  T)  is  the 
difference  between  the  ideal  reversible  double  layer  potential  for  concen¬ 
trations  R^,  and  the  actual  double  layer  potential,  therefore  it  repre¬ 
sents  the  sum.  effect  of  activation  and  concentration  polarization.  Even 
for  redox  reactions,  equation  2  is  not  the  only  form  of  rate  equation  which 
may  be  applicable  (see  4.5  later),  but  it  forms  a  good  starting  point  in 
the  lack  of  definite  knowledge  of  a  particular  reaction.  It  assumes  that 
the  solutions  considered  are  concentrated  enough  to  render  diffuse  double 
layer  effects  small  and  that  R  and  P  are  the  actual  reactant  and  product 
for  the  slow  step  in  the  reaction  (or  have  a  simple  equilibrium  relation 
with  the  actual  reactant  and  product) . 

3.2  Rates  of  Mass  Transport 

Let  us  consider  a  simple  anodic  redox  system  such  as 

j^Z+  electrode  p(Z+n)+  ^  ^3^ 

Associated  with  the  cations  will  be  the  corresponding  anions,  and  at  steady 
state,  the  law  of  neutrality  requires  that  each  step  in  equation  3  must  be 
accompanied  by  n  charge  transfer  in  the  electrolyte  between  the  electrodes. 
Each  cation  charge  in  the  electrolyte  must  be  balanced  by  a  corresponding 
anion  charge.  For  a  steady  state  condition  in  which  R  is  being  linearly 
transported  to  the  electrode  and  equimolar  amounts  of  P  are  being  transport¬ 
ed  away,  the  mass  transport  equations  for  a  stationary  system  can  be 
approximated  as  follows.  Let  D.be  diffusion  coefficient,  V  voltage  in  the 
electrolyte  and  x  distance  from  the  electrode  face.  Let  R  be  concentration 
of  reactant,  P  be  concentration  of  product  and  let  A  be  the  concentration 
of  anion  with  a  valence  charge  of  Z^.  If  the  mass  transport  effects  are 


such  that  different  parts  of  the  electrolyte  contain  different  totals  of 
R,  P  and  A  then  clearly  the  amount  of  water  present  per  unit  volume  is 
different.*  A  water  diffusion  gradient  exists  and  at  steady  state  this 
must  be  balanced  by  an  equal  but  opposite  flow  of  water^  that  is 

Djj^^dC/dx  -  Co  (4) 

C  is  the  concentration  of  water  and  o  is  the  bulk  flow  velocity  in  the 
opposite  direction  to  the  diffusional  effect.  This  bulk  flow  will  also 
carry  R,  P  and  A  with  it.  Whether  or  not  this  bulk  flow  significantly 
affects  the  transport  of  the  other  constituents^  R  for  example^  depends 
on  the  relative  magnitudes  of  D  dR/dx  and  Ro.  The  magnitude  of  D  dR/dx 
will  be  roughly  equal  to  that  of  D„  -dC/dx^  therefore,  from  equation  4 
the  effect  depends  on  the  relative  magnitudes  of  R  and  C.  It  will  be 
assumed  that  R^and  P,  A^  are  so  small  compared  to  C  that  the  bulk  flow 
of  water  carries  only  a  small  fraction  of  the  total  mass  transport  of 
the  constituents other  than  water.  This  will  be  a  fairly  good  approxi¬ 
mation  even  up  to  concentration  changes  of  several  moles  per  litre. 

It  will  also  be  assumed,  since  the  redox  ions  R  and  P  are  identical 
except  for  the  different  charge,  that  and  Dp  are  equal  and  the  ionic 
mobilities  (velocity  of  migration  In  unit  potential  gradient)  are 
proportional  to  the  valence  charge.  Then,  at  steady  states  rate  of 
transport  of  R  towards  electrode  (gm.  moles  per  sq.  cm.  per  sec.) 

X'j  ; 


*The  relation  Is  +  Xp  +  +  Xjj  q  ■  1  litre,  where  the  X  values  are  the 

effective  partial  volumes  of  the  constituents.  In  litres.  These  volumes 
are  functions  of  the  concentrations.  For  example,  the  value  of  X^  ■ 
where  Is  the  partial  molar  volume  of  R  (15)  In  litres  per  gm.  mole  per 
litre,  R  being  In  gm.  moles  per  litre.  The  functional  relation  of  to  R, 
P,  A,  C  Is  complex  and  can  vary  from  the  molar  volusie  of  the  pure  state 
to  zero  or  even  negative  values. 
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=  D  IdR/dx  +  RZu  dV/dx 

where  u  Is  the  mobility  per  unit  charge.  For  the  anodic  case,  dR/dx  will 
be  positive  and  dV/dx  negative.  Thus, 

1/nF  ■  D  dR/dx  +  RZiP  dV/dx  (5) 

Also,  considering  the  mass  transport  of  product  away  from  the  electrode, 

1/nF  “  -(D  dP/dx  +  P(Z+n)u  dV/dx)  (6) 

At  steady  state  the  total  transfer  of  the  non-dls charging  anion  Is  zero  and 

0  -  dA/dx  +  AZ^u^  dV/dx  (7) 

^A'  ^A  diffusion  coefficient  and  Ionic  mobility  per  unit  charge 

respectively  of  A.  algebraic  charge,  e.g.,  for  Cl^,  Z^  would  be 

-1.  The  final  equation,  that  of  electrical  neutrality  at  any  point  In  the 
electrolyte  Is 


ZR  +  (Z+n)P  -  A  (8) 

From  equation  8, 

ZdR/dx  +  (Z+n)dP/dx  •  Z^  dA/dx  (8a) 

Subtracting  equation  6  from  equation  5, 

(D/u)d(R+P)/dx  +  (RZ  +  P(Z+n))dV/dx  -  0 


From  equation  8, 


(D/u)  d(R+P)/dx  +  Z^AdV/dx  ■  0 


From  equation  7, 


(D/y)d(R+P)/dx  -  (D^/tJ^)dA/dx  -  0 
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Substltutlng  for  dA/dx  from  equation  8a, 

(D/M,)(dR  +  dP)-(D^/U,p[(Z/Z^)dR  +  ((Z+n)/Z^)dP]  =  0 
[  (D/u) -(W)  (Z/Z^)  ]dR  +  [  (D/ti)  -(D^/M^)  (Z4ii)  /Z^]dp:  -  0 
or 

dR  +  k2  dP  -  0  (9) 

The  boundary  conditions  are  R  ■  R^,  the  bulk  concentration  (unaffected  by 
mass  transport)  when  P  >  Pj^,  therefore 

k^  (R-Rj^)  +  k2  (P-Pi,)  -  0  (10) 

Equation  10  represents  the  relation  between  P  and  R  at  any  point  In  the 
electrolyte.  Eliminating  dV/dx  between  5  and  6, 

(l/nF  -  D  dR/dx)(l/RZ)  -  (-l/nF-D  dP/dx) (l/P(Z+n)) 

Substituting  for  dP/dx  and  P  from  equations  9  and  10  gives 

(1/nF-D  dR/dx)(l/RZ)  -  (-i/oF  +  D(kj/k2)dR/dx) (l/[ (k^/k2) (Ri,-R)  +  P^JCZ+n)) 

Collecting  terms  In  1 

1  -  (nFD  dR/dx)  (k2  +  k^k3)/(k2  +  k2k3)  (11) 

where 

k^  -  RZ/[(k^/k2)(R,j-R)  +  Pjj](Z+n)  (11a) 

Clearly  the  complete  form  of  equation  11  Is  too  complex  to  be  of  much  use. 

In  general,  however,  k^  will  be  some  value  between  0  and  oo  and 

1  -  (nFD  dR/dx) (1  to  k^/k2)  (12) 

k^/k2-  [(D/u)  -  (D^/y.^)(Z/Z^)]/[(D/u)  -  (D^/u^)(Z+n)/Z^J 


Now 
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But  (IB) 


D/|u  -  kT/e  “ 


(13) 


Therefore, 

^1/^2  =  (Z^-Z)/(Za-(Z+ii))  (14) 

Note  that  Z.  is  always  of  opposite  sign  to  Z  and  therefore  zero  cannot 

"  2+3+ 

occur  in  the  denominator  or  numerator.  For  example,  for  Fe  -»Fe  ,  with 

SO^  anion,  ;kj^/k2  =  (■4)/(-5)  =  0.8.  For  Sn^"*^ -»  Sn^^  with  Cl' anion, 

kj^/k2  “  (-3)/(-5)  =  0.6.  Thus  widely  different  from  1.  The 

above  equations  have  been  derived  assuming  that  Flck's  law  applies  with  a 

constant  diffusion  coefficient  and  that  mobilities  are  constant.  These 

assumptions  are  not  likely  to  be  exactly  true,  therefore  It  appears  that 

a  reasonable  simplification  of  equation  11  Is 


1  ■  nFD  dR/dx 


(15) 


where  U  Is  nearly  constant  with  varying  1  and  Is  of  the  same  order  of 
magnitude  as  the  true  diffusion  coefficient.  Equation  15  may  be  compared 
with  the  often  used  (17)  equation 

1  ■  (nFD/(l“n_|^))  dR/dx  (15a) 

n^  being  the  transport  number  of  the  reactant.  It  can  readily  be  shown 
that  this  equation  follows  from  equation  11  when  k^  Is  large  (Z+n  0,  for 

example),  Z  ■  and  D  ■  D^.  However,  If  equation  13  Is  taken  to  be 

precisely  correct,  D  ■  means  that  u  ■  u^  and  therefore  the  transfer 

number  Is  one-half.  Equation  15a  must  be  an  approximation  since  It  has  no 
term  In  the  diffusion  coefficient  of  the  anion. 

When  a  supporting  non-dls charging  electrolyte  Is  present  a  precise 
treatment  requires  the  writing  of  additional  equations  of  the  form  of 
equation  7  and  the  simultaneous  solution  of  the  equations  as  before.  How¬ 
ever,  It  Is  readily  seen  that  a  supporting  electrolyte  will  have  the  effect 


4 
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of  reducing  k^  towards  zero,  giving  i  =  nFD  dR/dx. 

For  diffusion  in  a  porous  medium  of  porosity  e  and  tortuosity.'  factor 
q  (18),  the  mass  transport  based  on  the  total  area  is  given  by 

i  =  nFD  --dR/dx 
err 

(16) 

where 

Deff  -  Ce/q 

(16a) 

3.3 

Redox  Reactions  at  a  Porous  Electrode: 

Ohmic  Effects  Neslected 

3.'3.1  Reaction  Within  Electrode 

The  system  considered  is  illustrated  in  Figure  1.  It  will 
be  assumed  that  electrolytic  current  and  mass  transfer  occur  through  one 
plane  face,  with  negligible  edge  effects  and  negligible  non-uniformity  of 
current  density  in  the  electrolyte.  The  system  is,  therefore,  one  dimensional 
and  we  can  consider  a  unit  area  of  face.  Let  the  concentration  of  reactant 
and  product  in  the  bulk  of  the  electrolyte  be  R^,  respectively.  Let  the 
thickness  of  the  electrode  be  t.  Equation  2  can  be  written  as 

or 

1  -  i^((R/R^)e°^'''’  -  (P/P^J)e■^^■^^^^^  (17) 

For  electrodes  of  reasonable  thickness  and  pore  size,  the  macro-pore  size 
will  be  at  least  two  orders  of  magnitude  less  than  the  thickness  (10  microns 
and  1  mm,  for  example)  and  we  can  assume  that  mass  transfer  effects  across 
the  radius  of  a  pore  are  negligible  compared  to  mass  transport  effects  along  the 
axis-  of  the  pore. 

Consider  e  differential  element  dx  as  shown  in  Figure  1.  The  reacting 
area  is  Sdx  and  the  current,  dl,  from  the  element  is 

dl  -  io((R/R^)e°^^^  -  (P/Pj^)e"^^‘“^^^'’)Sdx 
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The  amount  of  reactant  diffused  out  of  the  element  at  the  plane  x  dsID'dr/dx 
Xwhere  D  is  understood  to  be  the  of  equation  16).  The  amount  diffused 

in  through  the  plane  at  x  +  dx  is  equal  to  D(dR/dx  +  d(dR/dx)).  In  the  usual 
way^  therefore,  at  steady  state, 

di  ■  nFD(d\/dx^)dx 

and 

d^R/dx^  =  (i^S/nFD)[(R/R^)e°^''^  -  (18) 

If  ohmic  effects  are  neglected,  T]  is  constant  through  the  electrode.  For  a 
redox  reaction  with  supporting  electrolyte  the  effective  diffusion  coefficients 
of  R  and  P  are  nearly  equal  and  little  error  is  Introduced  by  letting  the 
diffusion  coefficients  be  equal.  Then,  for  equl -molar  counter. -ciitrent  dif¬ 
fusion 

R  +  P  -  +  P^^  (19) 

Substituting  for  P  into  equation  18 

d^/dx^  -  (l^S/nFD)[(e“^^^/R^)  +  (e"^^"“^’^^'’/P^)  jR 

-(l^S/nPD)(R^j  +  P^)(e‘^^'^^^^*’/Pjj)  -  K^R  -  K2  (20) 

The  boundary  conditions  are  dR/dx  ■  0  at  x  ■  0,  since  no  mass  transport  occurs 
beyond  the  left  hand  face,  and  R  ■  R  at  x  -  t,  R  being  the  concentration  at 
the  right  hand  face.  The  solution  of  equation  20  for  these  boundary  conditions 

is 

K1R-K2  _  e'^*  +  e 

I  .e^t 

The  total  steady  state  current  from  reaction  in  the  Interior  of  the  electrode 

is 
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i  -  nFD(aR/dx)^_j. 

Now 

K,®/dx - -  (  ^/l^  e  'Si*  -  ^  e"  ■Si*) 

therefore 

(K,R  (.'Si'  -e-'STt) 

i  .  nFD - - - 

(Kp  (  ^ 

=  HFOCK^Rg-KpClZ-yK^)  tanh(ts/^) 

But  since  the  RHS  of  equation  20  equals  the  RHS  of  equation  18^  Kj^R^-K2  can 
be  replaced  and 

1  =  Cl^St)[(R^/R^)e°^/'’  -  (Pg/P^^)e'^^"^^^/’’](l/t  J^)tanh(t*/K]^)  (21) 

Thus  we  arrive  at  the  usual  form  of  rate  equation  but  modified  by  a  factor 
(  1/Wk^)  tanh(  t  »/k^)  . 

Several  limiting  cases  can  be  described. 

(1)  Low  exchange  current,  low  to  moderate  polarization.  From  equation 

20, 

=  (i^S/nFD)[(e°^^^/Rj^)  +  ]  (22) 

If  l^S  is  low  and  polarization  not  too  high,  t  >Jk^  Is  small.  Then 
(1/t  VKj^) tanh  t  *yKj^  is  near  1,  and' 

1  -  (i^St)[(Rg/R^)e°^/‘’  -  (P^/P^)e'^^'®^^/*’]  (23) 

Thus,  If  there  Is  negligible  concentration  di^p  from  the  bulk  electrolyte  to 
the  electrode  surface  (as  would  be  expected/at  low  currents),  R^  ■  Rj^,  P^  ■  Pj^ 


*For  example,  (tanhx)/(x)  ■  0.97  at  x  >  0.3. 
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and  we  obtain  the  usual  form  for  a  low  exchange  current  reaction  with  the 
whole  of  the  interior  surface  area  (St)  of  the  electrode  being  utilized. 

If  ipSt  is  so  small  that  t  »>/k^  is  small  for  polarizations  above  50  milli¬ 
volts  or  so,  then  equation  23  reduces  to  the  normal  Tafel  form,  with  a  slope 
of  2.3  RT/CtoF. 

(ii)  Low  or  moderate  exchange  current,  high  polarization. 

From  equation  22  it  is  seen  that  as  T]  goes  up  (or  as  i^S  is 
higher),  t  goes  up.  Eventually  the  situation  is  reached  where 
t n/k^  >  2.  Then  tanh(t «/k^)  is  equal  to  one  and  equation  21  goes  to 

i  -  i^S  [(Rg/R,j) 

N/(i^S/nFD)[(e°^^^/R^)  -  e'^^"“^'^^‘’/Pj^)  ]  (24) 

The  thickness  of  the  electrode  disappears  from  the  equation  since  the 
concentration  of  reactant  in  the  interior  reaches  zero  at  some  penetration 
less  than  the  thickness.  When  T]  is  great  enough  for  e  to  be 

negligible, 

i  =  (i^SnFDR^)  ^  (Rg/Rj^)e°^/^^  (25) 

For  negligible  concentration  gradient  outside  of  the  electrode  R^/R^,  =  1 
and  equation  25  is  a  Tafel  form  equation,  with  a  slope,  however,  of  twice 
the  expected  normal  value.  The  apparent  exchange  current  is  »/T^rJDR^. 

(iii)  High  exchange  current;  pure  concentration  polarization. 

For  very  high  exchange  currents,  that  is,  i/i^St  small,  we  get 
the  normal  form 

T]  -  2.3b  log  (R^^Pg/RgP^)  (26) 

The  physical  meaning  here  is  that  the  current  can  be  supplied  from  an 
infinitesimal  layer  at  the  electrode  surface.  This  layer  is  so  thin  that 
the  concentration  gradient  which  is  diffusing  reactant  from  the  surface 
into  this  reaction  zone  is  virtually  zero.  The  interior  of  the  electrode 
is  at  Rg'SR  and  +  &F;  where  5R  and  &F  are  small  quantities  which  bring 
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the  interior  of  the  electrode  to  equilibrium  conditions.  In  the  infinitesimal 

surface  reaction  layer  where  ?.  goes  from  R.-5R  to  R  ,  the  slight  deviation 

s  s 

from  equilibrium  is  sufficient  to  supply  the  current  being  drawn. 

3.3.2  Reaction  at  Surface  of  Electrode 

The  treatment  in  the  previous  section  describes  the  current - 
voltage  relation  for  current  produced  from  the  interior  of  the  electrode.  In 
addition,  however,  the  exposed  face  of  the  electrode  will  support  current. 
Although  the  geometric  external  area  considered  is  unit  area,  the  external 
reacting  area  which  is  freely  available  without  mass  transport  restrictions 
in  the  interior  porous  system  can  be  considered  as  £  sq.  cm.  per  sq.  cm.  of 
geometric  area  (see  discussion  of  equation  1) .  Therefore,  a  more  complete 
equation  for  current  density  is 


^  "  ^interior  ^exterior 


-  (l^tS)[(Rg/Rjj)e°^/‘'  -  (P^/P^^)e■^^■°‘^^/''](l/t^/K[)tanh(t^/K5^) 

+  i^f  (!-€)[  (R^/R^)e“^^'’  - 

l-€  is  the  area  of  the  exterior  allowing  for  open  pore  mouths.  Thus, 

1  -  [(i^tS)(l/t*/Kj^)(tanht»/I^)  +  l^f  (1“G) ] [ (R^/R^)e°^^’’  -  (Pg/P^,)e■^^‘®^^^'’J 

(27) 


The  relative  importance  of  the  two  currents  is  described  by 


h„terlor/‘.»t.rIor  '  V5^)t.nh(t 

From  equation  1, 

W^^ext  “  [(t(l-€)A](l/t«;^)tanh(t^)  (28) 

A  is  the  Internal  perimeter  area.  For  example,  an  electrode  of  30%  porosity 

3 

and  a  pore  diameter  of  10  microns  has  A  ■  4  e/d  ■  (4) (0.3) 10  %  1000  sq.  cm. 
per  cubic  cm.  For  t  small  (<  0.3), 
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For  a  thickness  of  0.1  cms.  and  A  ■  1000  cm.  the  ratio  of  interior  to 
exterior  current  is  then  about  100^  and  the  exterior  current  is  negligible. 
For  t  >  2, 

h«'‘ext 

1/2 

-  A(l-€)/*y(l^S/nFD)  ] 

(28b) 

For  negligible 

From  equation  (25) ^ 
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.  1  +  7i/l^ 

As  before  It  Is  assumed  that  the  effective  mass  transport  coefficients  of 
R  and  P  are  equal.  7  is  the  ratio  R^j/P|j-  Substituting  in  equation  21, 

i  -  (i^t8)(l/^/K]^)tanh(t^/i^  [(l-l/ij^)e“^^^  -  (1  +  7l/ij^)e’^^““^^^^] 

Rearranging; 

(Its/.)  (.»l/b  .  ^-(l-a)V>>) 

i/i'  m  O  h  _  _  . 

^  (t^/K^/tanh(t^/K^)  +  (i^tS/i^)(e°^^’^  +  (30) 

A  somewhat  more  convenient  form  can  be  obtained  as  follows.  From  equation  22 

-  (l^tS/tnFDR^)[e“'^^'’  + 

The  limiting  current  is  given  by 

ij^  =  nFDRj^q/6€  (31) 

8  is  the  effective  thickness  of  the  "stagnant  film"  at  the  exterior  of  the 
electrode,  q  and  €  are  present  to  convert  the  effective  diffusion  coefficient 
in  the  porous  electrode  to  that  in  the  free  electrolyte  outside.  Let 


tl/  -  at/de  ■  ease  of  mass  transport  at  exterior 
ease  of  mass  transport  in  interior 


Also  let  s  ^o^^o  exchange  current  density  for  complete  utili¬ 

zation  of  the  electrode.)  Then; 


Equation  30  now  goes  to 
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i/ij  - 


(i  /i,)(e«n/b  . 

O  L 


ta„h  '‘fd  (33) 

O  L 

Equation  33  represents  the  complete  form  up  to  the  limiting  current.  When 
1  is  not  near  1^^  It  reduced  to  the  previous  forms,  with  Rg/R^  “!• 

For  high  exchange  current  reactions,  that  is,  large,  -♦  0> 

equation  33  reduces  to 

i/1,  -  -  e-(l’«>n/b)/(gOflr)/b  ^  ^^-(l-a)Tl/bj 

L 


Rearranging, 


Ti  =  2.3b  log  [(1  +  7i/i^)/(l-i/i.)] 


This  is  the  expected  equation  for  pure  concentration  polarization. 

3.3.4  Calculation  of  Equations 

Equation  33  contains  the  parameters  i  /i  ,  ip,  y,  a,  b.  For 
irreversible  conditions,  y  disappears  and  a,  b  can  be  combined  as  Qt/b.  It 
will  be  noted  that  for  irreversible  conditions  (i^/i^^)  and  always 

occur  together  as  Thus  as  (i^/i^^)  varies,  the  shape  of  the 

curve  log  (i/ij,)  versus  CWl/b  will  be  the  Same,  but  displaced  by  a  definite 
Increment  of  arj/b.  If  we  know  the  shape  at  (i^/ij^)]^  and  if  At)  represents 
the  bodily  shift  at 

(I  /I  =  (i  /i,),  e<"1l/l>)  +  (ctl-l/l.) 

'  o  L  1  '  o  L  2 

or  ATI  =  (2.3)(b/a)  log[(i^/ij^)^/(i^/ij^)2]  (35) 

It  is  convenient  to  compute  the  results  in  terms  of  i/1^  and  crn/b  ■  T)  . 

Then, 


-(l-a)'n/b  -  Tj[(l-a)-l] 
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Reasonable  values  of  ^  can  be  derived,  as  follows.  6  is  normally  about  0.05 

cms  for  unstirred  electrolytes  and  may  be  10  times  lower  for  vigorous  rates 

of  flow  past  the  electrode  (17))..  q/e  is  normally  in  the  range  5  to  10,  t 

may  be  from  about  0.01  cms  (a  lower  value  may  mean  that  the  external  surface 

reaction  cannot  be  neglected) tto  say,  1  cm.  Therefore 

qt/6€  -  (5)  (0.01) /( .05)  ■  1  minimum 

or  »  10/ (.005)  =  2000 

For  unstirred  electrolytes,  q/€  =  5  and  a  thickness  of  0.1  cm.,  a  value  of 
^  s  10  seems  reasonable.  Therefore  results  were  computed  for  ^  1,  10,  100, 

For  near  reversible  conditions,  y  was  taken  as  1,  OC  as  1/2.  The  results  are 
shown  in  Figure  2. 

[Note  that  f  =  0  means  t  =  0.  When  t  =  0,  the  term  for  reaction  at  the 
external  surface  must  be  included  and  current  from  the  interior  neglected. 
Therefore  tp  =  0  represents  a  solid  electrode  with  a  surface  exchange  current 
density  of  i^.]  It  can  be  seen  from  Figure  2D  that  as,. jjj  becomes  greater,  a :Tafel 
section  is  obtained  with  double  the  normal  slope  and  this  section  increases  in 
length  as  tp  becomes  greater.  However,  it  should  be  noted  that  when  the 
reaction  is  Irreversible,  not  near  the  limiting  current  and  completely  in  the 
Internal  mass  transport  affected  region  ( */^  large),  equation  33  (or  equation 
25)  reduces  to 

Since  and  Tp  are  both  proportional  to  t,  it  is  not  possible  to  force  the 
system  into  a  longer  Tafel  region  by  increasing  t.  If  the  basic  exchange 
current  is  so  high  that  no  clear  Tafel  region  exists,  then  increasing  t  will 
not  produce  a  Tafel  region.  The  apparent  exchange  current  of  this  region  is 
i.  (i  /i.  )  /Tp.  For  the  case  plotted  in  Figure  2D,  when  ^  ■  100  this  equals 

Lf  O  Li 

1^/100,  which  agrees  with  the  value  obtained  by  extrapolation  of  the  linear 
part  of  the  curve  to  zero  polarization. 
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3.4,  Redox  Reactions  at  a  Porous  Electrode;  Ohmic  Effects  Included. 

3.4.1  Basic  Equations. 

The  treatment  given  in  section  3.3  is  only  a  limiting  case 
of  the  geherat  current -polarization  relation,  since  it  neglects  the  effect  of 
ohmic  voltage  gradient  in  the  pore  electrolyte.  The  effect  of  potential 
gradient  in  the  electrolyte  within  the  pores  of  the  electrode  depends  on  the 
relative  magnitudes  of  the  conductivity  of  the  electrolyte  and  the  diffusion 
coefficient  and  concentration  of  the  reactant.  For  a  high  concentration  of 
supporting  electrolyte  (and  consequently  high  conductivity)  and  a  low  concen¬ 
tration  of  reactant,  the  limiting  current  will  be  low  and  the  polarization  is 
controlled  mainly  by  activation  and  mass  transport.  For  a  weakly  conducting 
electrolyte  and  a  high  concentration  of  reactant,  the  activation  and  the 
ohmic  effects  within  the  pores  will  be  the  controlling  factors. 

Let  us  consider  the  case  where  the  conductivity  of  the  electrolyte  is 
reasonably  constant,  even  though  reactant  ions  are  being  changed  into  product 
ions.  This  is  a  good  assumption  for  a  strongly  supported  electrolyte.  Let  p 
be  the  specific  resistance  of  the  electrolyte  in  the  pores.  Then,  as  before, 

p  =  p^€/q  (36) 

where  p^  is  the  free  specific  resistance.  Ohmt  law  gives 

i  =  (l/p)dV/dx  (37) 

When  current  is  being  drawn  from  the  electrode,  the  current  is  supported  by 
ionic  flow  in  the  electrolyte  in  the  pores.  Just  as  in  the  case  of  mass 
transport,  all  of  the  current-carrying  ions  must  pass  through  the  external 
face, while  only  a  few  may  be  required  to  reach  the  left  hand  face.  The 
voltage  gradient  is  highest,  therefore,  at  the  external  face.  Unlike  ohmic 
loss  in  the  electrolyte  away  from  the  electrode,  the  voltage  change  in  the 
pore  electrolyte  has  a  direct  effect  on  the  speed  of  the  electrochemical 
reaction,  since  the  pore  electrolyte  potential  is  the  potential  from  which 
ions  are  discharging  to  the  surface.  In  the  equation 
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the  polarization  must  be  measured  between  the  pore  electrolyte  and  the  electrode, 
T]  will  vary  along  the  pore,  being  greatest  at  the  external  face,  and  dV/dx  ■ 
dr\/dx. 

The  equations  to  be  solved  are  now 

di  »  i  S((R/R.  -  (P/P,  )e’^^'°‘^''^^^)dx 

O  □  D 

dl  “  nFD(d^R/dx^)dx 
i  «  (1/p)  dTl/dx  , 

The  boundary  conditions  are 

dR/dx  =  0  at  X  =  0 

R  =  R  at  X  =  t 
s 

T]  =  T)s  at  X  =  t 
T)  =  at  X  =  0. 

The  total  interior  current  is  given  by 

i  =  D(dR/dx)^^^  =  (l/p)(dTi/dx)^j.. 

As  before,  we  can  assume  the  mass  transport  factors  of  R  and  P  to  be  equal  and 
constant  and 

R  +  P  =  Rb  + 

or  P/Pj^  =  (7  + 

Combining  the  ohmic  and  mass  transport  equations 

(l/p)dTl/dx  ■  nFD/(d^/dx^)dx 
o 

Since  dR/dx  ■  0  at  x  -  0, 
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fdr\  "  ptiFD  /dR 

T]  -T]  -  pnFD  (R  • -R)  (38) 

s  s 

We  are  now  in  a  position  to  compare  the  relative  effects  of  mass  transport 

and  internal  ohmic  voltage  gradient.  For  example,  if  R  =  0  at  x  =  0, 

“3  3  •62 

R  =  10  gm.  mole/cm.  ,  p  =  50  ohm  cm.,  n  =  1,  D  =  10  cm.  /sec.,  then 
s 

^s'%  2;  millivolts.  Such  a  small  change  in  voltage  will  not  produce 

much  effect  compared  to  the  corresponding  depletion  of  reactant  in  the 

interior  of  the  electrode.  As  current  density  is  increased,  the  penetration 

into  the  electrode  at  which  R  becomes  nearly  zero  becomes  smaller  and  due 

to  this  mass  transport  effect  the  reaction  is  concentrated  in  the  region 

towards  the  right  hand  face  of  the  electrode.  Then  the  mean  path  through 

which  ions  have  to  travel  becomes  less  and  although  i  becomes  greater  and 

hence  dT]/dx  becomes  greater,  this  is  accomplished  by  x  being  smaller  and 

not  by  T]  becoming  greater;  'ng“'n  remains  the  same. 

T]  -T)  becomes  greater  as  p  and  R  become  greater.  However,  if  R  is 
s  s  s 

an  ion  it  is  not  possible  to  make  the  specific  resistance  p  greater  while 

increasing  or  maintaining  R  .  To  make  the  internal  voltage  effect  over- 

s 

weigh  the  mass  transport  effect  it  would  appear  to  be  necessary  to  use  low 
conductivity  solutions  with  a  high  concentration  of  a  non-conducting 
reactant. 

The  complete  equations  can  be  written  as  follows.  From  equation  38 
R  -  Rg  -(T]g-Tl)/pnFD  (39) 

Combining  the  kinetic  and  ohmic  equations  gives 

(l/p)d^Tl/dx^  -  i^S((R/R^)e°^^^  -  (P/P^)e'^^‘°'^^^'’) 

-i^S(y  + 
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Substituting  for  R  from  equation  39 

(l/pi^S)d^1l/dx^  -  [(e°^^''/Rjj)  +  (e‘^^'®^^^^/Pj^][Rg-(Tig-Tl)/piiFD] 

-  (7  +  (40) 

Equation  40  Is  the  equation  relating  polarization  to  distance  through  the 
electrode. 

3.4.2  Limiting  Cases. 

Case  1.  Low  resistance,  low  current  density,  R  small. 

This  would  apply  for  a  low  concentration  of  reactant  and  a  high 
concentration  of  supporting  electrolyte.  Under  these  circumstances,  the  mass 
transport  limited  currents  would  be  low;  the  limit  at  which  ohmic  effects  are 
negligible  is  the  case  which  has  been  dealt  with  in  section  3.3.  Mathematical¬ 
ly,  from  equation  38  this  limit  is  given  by  being  small,  for  low  values 

of  p  and  Rg.  This  case  is  probably  the  most  widely  applicable. 

Case  2.  High  exchange  current:  pure  concentration  polarization. 

For  very  high  exchange  currents,  a  slight  deviation  from  the 
reversible  potential  will  give  a  large  current.  Therefore,  as  discussed  in 
section  3.3.1,  the  reaction  will  be  accomplished  in  a  thin  layer  at  the  out¬ 
side  of  the  electrode.  As  this  layer  tends  to  zero  thickness,  that  is 
dx  0,  both  dr  and  dr|  tend  to  zero  and  the  electrolyte  in  the  interior  of 
the  electrode  is  at  almost  the  same  concentration  and  potential  as  at  the 
surface.  Equations  26  and  34  then  apply  at  this  limit. 

Case  3.  Very  low  conductivity,  high  reactant  concentration,  high 

current. 

This  case  applies  when  the  currents  being  drawn  are  great  enough  compared 
to  the  conductivity  and  small  enough  compared  to  the  mass  transport  that  ohmic 
polarization  is  large  without  the  values  of  R  being  much  changed.  Mathematical¬ 
ly,  from  equation  38,  -T]  must  be  large  for  R  -R  small.  For  example,  a  change 

of  R  from  Rg  to  0«9Rg  will  not  produce  much  change  in  rate  or  polarization. 

The  change  in  polarization  through  the  electrode  is  then 
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T|  "T|  ^  pnFD(R  /lO) 
s  s 

*"3  3  "6  2  "1 

For  p  =  10000  ohm  cm.,  R  =  2N  =  2(10  )  gm.  moles/cm.  ,  D  «=  10  cm.  sec.  , 

s 

T)  -Ti  would  be  0.2  volts.  Such  a  voltage  change  would  completely  overweigh.: 

8 

the  effect  of  change  in  R.  The  equation  to  be  solved  is  then 

d^Tl/dx^  -  pi^S[(R^/Rj^)e°^/^  -  (41) 

The  physical  picture  Is  as  follows.  At  the  left  hand  face  of  the  electrode, 
with  a  polarization  r|^,  the  current  being  produced  is  very  small.  However, 
this  small  current  produces  a  small  ohmic  gradient  which  Increases  T)  as  we 
precede  from  left  to  right.  The  larger  produces  a  larger  current  increment 
and,  since  the  total  current  is  additive,  the  ohmic  effect  becomes  greater 
and  greater.  The  exponential  nature  of  the  current  increment  versus  polari¬ 
zation  means  that  at  some  cumulative  polarization  the  current  will  increase 
extremely  rapidly  and  most  of  the  reaction  occurs  towards  the  right  hand  face 
of  the  electrode.  It  is  clear  that  very  small  changes  in  will  produce  very 
large  changes  in  r|^  and  the  final  current.  Thus  we  may  cover  several  orders 

of  magnitude  of  current  with  an  almost  constant  value  of  n  .  The  value  of  Ti 

*o  *o 

for  a  given  order  of  magnitude  of  current  will  depend  on  Si^.  Again,  two 
limiting  cases  can  be  considered. 

Case  3A.  Low  exchange  current. 

Let  us  consider  the  case  where  the  effective  exchange  current  Si^ 
is  so  low  that  must  be  fairly  large  before  appreciable  current  is  obtained. 
For  this  irreversible  condition,  the  negative  exponential  term  in  equation 
41  can  be  neglected  and 

d^Tl/dx^  -  pi^S  (Rg/Rj,)  (42) 

Letting  'H  “  1)^  + 

d^(A'ri)/dx^  =  pi  S(R  /R.)e“^o^'’  (42a) 

O  S  D 

The  boundary  conditions  are  now  At^  ■  0  at  x  ■  0,  d(AT))/dx  -  0  at  x  ■  0.  For 
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very  low  values  of  pi  S(R  /R,  the  current  is  small,  ohmic  effects  are 

2  s  b  ^ 

negligible  and  d  AT]/dx^  -♦  0.  Equation  42a  can  be  put  as 

d(ATl/dx)  =  pi^S(Rg/R^)e“^o/bj-'^eQ!^'n/bjjjj  . 

'■‘•v  ^ 

2  ‘2  ■  ■ 

When  d  A^/dx  -»  0,  AT)  -»  0, 

=  pi^S(Rg/R^)e°^o/‘’x 
Since  i  *  (1/p)  (dT)/dx) 

X®C  f 

i  =  i  St(R  /Rje°^o/'’ 
o  '  s  b' 

This,  of  course,  is  the  required  Tafel  form  with  uniform  utilization  of  the 
electrode.  The  more  interesting  case  is  when  pi^S  is  high  enough  to 

make  At]  appreciable,  so  that  the  bulk  of  the  reaction  occurs  in  a  layer  near 
the  external  face.  Putting  pioS(Rg/R^)e*^o^'’  ■  k  and  a/b  ■  a,  the  equation 

d^(ATl)/dx^  -  ke®^'^ 

Let  y  ■  d(Al))ydx  then  d^Aq/dx^  ■  dy/dx  ■^y/dAT))y. 

Then 

ydy/d(AT])  =  ke^^ 

and  y^/2  =  (k/a)e^^  +  C 

or  dAT]/dx  *  [(2k/a)e^'^  + 

When  X  «  0,  At]  =  0  and  d(AT))/dx  =  0,  therefore  C  =  -(2k/a)  and 

d(ATi)/dx  =  [(2k/a)e®^  -  (2k/a)3^^^  (43) 

Thus  the  total  current  is 

i  -  (l/P)[d(/!in)/dx]^„,. 

=  (1/p)  ^  (2k/a)  (e®^'^s-Tlo)  _  i)l/2 

.  (l/p)^/2pi  S(R  /R  )  e‘^°^'’b/a  (g^^’^s/'Ho) /b  _  ^^1/2 
o  so 


(44) 
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When  exp[a(T)  -Ti  )/b]  is  large  enough  for  the  1  to  be  negligible^ 

8  O 

/ 

.  ..  /  2i  S(R  /Rjb  ^ 

i  a  /  °  ^ ^ (45) 

V  a  p  ^ 

For  Rg/Rjj  ■  the  plot  of  log  i  versus  T)^  is  a  Tafel  form  with  double  the 
usual  slope.  The  apparent  "exchange  current"  is  dependent  on  p. 

Case  3B.  High  exchange  current. 

Letting  “  1^  ^s^^b  "  equation  41  can  be  written  as 

d^T]/dx^  =  sinh(ar)/b)  (46) 

This  has  been  solved  by  Stender  and  Ksenzhak  (3).  (See  also  Section  4.2). 

For  the  boundary  condition  of  0^  (that  is^  high  exchange  current  and 

thick  electrodes) 


i 


/ 


^^pSb  .an /2b 

ap  '' 


g-(l-a)n/2b^ 


(47) 


We  will  not  discuss  equations  45  or  47  at  this  point,  because  the  con¬ 
ditions  under  which  they  apply  will  be  very  rarely  met  in  the  type  of  system 
considered  here.  In  this  respect  we  disagree  with  Stender  and  Ksenzhak,  who 
Calculated  results  at  p  =  50  ohm-cm,  assuming  uniform  concentration  of  reactant 
throughout  the  electrode.  Equation  38  shows  that  would  be  only  a  fraction 

of  a  millivolt,  with  complete  internal  mass  transport  control  of  the  process 
for  such  a  resistance  and  a  feasible  reactant  concentration.  For  the  same 
reasons,  we  will  not  compute  equation  40  since  it  is  felt  that  the  ohmic  effect 
will  be  small. 

Another  interesting  case  is  that  where  no  supporting  electrolyte  is 
present  and  the  product  is  non-ionic.  The  specific  resistance  is  then  given 
by 


P  -  P.<R./lt) 

where  p  is  the  resistance  at  concentration  R  . 

S  8 


(48), 

For  1  N  solution,  R  ■  10  ^ 
*  s 
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3 

moles/cm  and  p  is  about  50  ohm  cm  maximum. (allowing  for  porosity,  tortUosity 
8 

factors) .  The  ohmic  and  mass  transport  equations  now  combine  to  give 

dT)/dx  =  p  nFDR  (l/R)/(d^R/dx^)dx 
s  s 

o 

Since  dR/dx  ■  0  at  x  =  0, 

/dT)  =  p  nFDR  /(l/R)dR 

8  S 

T)  -  Tl  -  p  nFDR  2.3  log(R  /R)  (49) 

S  8  8  8 

“2  *“6 
Since  p  R  is  about  5(10  ),  D  is  about  10  , 

8  8 

\  2;  (1/100)  log(Rg/R) 

10  log(R  /R)  millivolts, 
s 

A  value  of  R  /R  of  100  would  represent  large  mass  transport  effects  in 
8 

the  pore  electrolyte,  but  would  be  only  20  millivolts.  Therefore, 

although  the  ohmic  effect  is  much  greater  in  this  case,  it  may, not  cause  a 
very  pronounced  deviation  from  the  treatment  of  section  3.3..  It  shbiildsbe 
noted  that  at  a  given  current  density,  the  treatment  of  section  3.3  gives  a 
value  of  T)  which  lies  between  and  1)^.  Therefore,  the  deviation  of  t) 
from  is  some  fraction  of  T]^  -  not  the  total  amount.  The  complete 
equations  can  be  set  up  as  follows.  The  combined  mass  transport -electro¬ 
chemical  equation  is  equation  18 

d^/dx^  -  (i  S/nFD)[(R/R.  )e“^'^^  -  (P/P, 

O  D  D 

The  measured  polarization  is  that  at  the  surface,  r|  and  from  equation  49, 

s 

=  (R/R  )aPs"^DR^/b  aris/h 
s 

Substitution  of  this  into  the  preceding  equation  gives  the  final  form.  The 
effect  of  ohmic  voltage  gradient  is  given  by  the  factor 
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(R/r/P.'’™.'" 

For  a  ■  1/2,  b  “  0.026  at  room  temperature  and  p  NFDR  ~  5/1000  we  get 

S  8 

(R/R  ^  )1/10 

s  s 

At  R/R  =  10  the  factor  is  10  s  0.93,  therefore,  the  deviation 

s 

from  1  is  not  great  and  the  solution  of  equation  18  as  before  can  be  used. 

Although  T)  -T)  may  be  appreciable,  Tj  lies  close  to,  but  a  bit  greater  than, 

the  value  of  T]  calculated  from  equation  33.  The  physical  meaning  here  is 

that  the  ohmic  effect  is  small  until  R/R  is  small,  that  is,  well  within  the 

s 

electrode.  But  most  of  the  reaction  occurs  towards  the  outer  face  of  the 
electrode  so  that  ‘Hg  %  11  effective.  may  be  markedly  different  to  T)^,  but 
very  little  reaction  occurs  at  conditions  due  to  the  low  value  of  R. 

3.5  Effect  of  Double  Electrode  Structure 

The  previous  treatments  have  assumed  that  the  electrode  is  homo¬ 
geneous  with  penetration  and  that  the  polarization  is  measured  adjacent  to 
the  exposed  face  of  the  electrode.  If  a  double  layer  structure  is  employed, 
as  shown  in  Figure  3,  the  mass  transfer  coefficient  in  the  outer  layer  is 
different  to  that  in  the  free  electrolyte.  Let  us  suppose  that  the  inner 
layer  is  the  active  electrode  material,  while  the  outer  layer  is  inert 

electrochemically.  The  polarization,  and  the  value  of  R  ,  to  be  used  in 

s 

the  previous  equations  are  those  at  t,  whereas  the  polarization  is  measured 
at  L. 

If  €/q  is  the  porosity-tortuosity  factor  for  the  coarse  layer,  and  5  is 
the  effective  "film  thickness"  in  the  free  electrolyte, 

i  =  nFD(Rjj-Rg)  l/[(Lq/e)  +  6] 

=  nFDRj^/[(Lq/e)  +  6]  (50) 

and 

R^/R.  =  1-i/l^ 

S  D  Li 

This  is  the  result  used  previously,  except  that  1^  now  depends  on  the  values 
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of  L,  q  and  e  of  the  outer  layer. 

For  example,  if  L  is  0.05  cms  and  q/e  is  5  (a  reasonable  figure),  the 
term  Lq/€  is  0.25  cms.  This  is  five  times  the  effective  6  of  0.05  cm.  for 
unstirred  electrolytes,  therefore  i^^  will  be  only  1/6  of  the  usual  limiting 
current. 

The  ohmic  voltage  drop  in  the  coarse  layer  is 

An  =  ip'Lq/e  (51) 

The  measured  polarization  must  be  corrected  for  this  An  before  experimental 

results  are  compared  to  theoretical  expectations.  For  a  large  conductivity, 

2 

p'q/e  2;  10  ohm  cm  and  for  L  =  0.05  cm,  i  =  100  mamps/cm  ,  An  =  50  millivolts. 

It  is  apparent  that  the  use  of  a  double  layer  structure  is  a  big  dis¬ 
advantage  where  the  limiting  factor  is  the  mass  transport  of  reactant  from 
the  bulk  of  the  electrolyte  to  the  reacting  electrode. 

(4)  GAS  DIFFUSION  ELECTRODES  (THREE  PHASE  SYSTEMS) 

4.1  General  Introduction 

The  mode  of  operation  of  a  gas  diffusion  electrode  is,  in  general, 
more  difficult  to  describe  mathematically  than  a  simple  redox  system.  The 
principal  difficulty  is  that  the  processes  occurring  at  the  three  phase  boundary 
are  not  known  with  any  degree  of  certainty.  It  is  necessary  to  choose  a 
simplified  model  of  what  we  think  is  happening.  If  an  attempt  is  made  to 
allow  for  every  possible  factor,  the  model  becomes  so  complex  that  extensive 
computation  is  required.  It  would  appear  to  be  more  reasonable  to  start 
with  the  simplest  models,  determine  how  these  would  behave,  and  compare  the 
theory  to  experimental  values.  For  example,  a  simple  model  which  might  be 
applicable  under  certain  conditions  is  shown  in  Figure  4.  In  this  model  it 
is  assumed  that  the  bulk  of  the  electrolyte  is  stabilized  in  the  macro-pore 
system  of  the  electrode,  but  that  a  thin  film  of  electrolyte  covers,  and 
penetrates  into,  the  micropore  system.  Reaction  occurs  by  diffusion  of  gas 
through  the  film^followed  by  adsorption  and  reaction  on  the  active  surface. 

The  current  must  be  supported  by  movement  of  ions  through  the  thin  film  of 
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electrolyte,  parallel  to  the  active  surface.  If  the  "film"  is  so  thick  that 
it  fills  the  pore,  then  gas  will  not  be  able  to  diffuse  to  the  active  surface 
fast  enough  to  support  much  current.  On  the  other  hand,  if  the  film  is  very 
thin,  the  horizontal  ionic  mass  transport  and  conductivity  in  the  film  will 
be  too  low  to  support  much  current.  Thus  this  model  can  only  apply  if  the 
film  thickness  is  considerably  less  than  the  macropore  size  (10  microns,  say) , 
but  much  thicker  than  a  few  molecules. 

The  treatment  of  such  systems  may  sometimes  reduce  to  similar  forms  to 
those  considered  before,  but  it  might  be  expected  that  fresh  formulae  will 
have  to  be  derived  in  some  cases. 

4.2.  Thin  Film  Model;  Hydrogen  Fuel  in  Acid  Electrolyte. 

The  case  of  hydrogen  fuel  used  at  a  porous  diffusion  electrode 
in  acid  electrolyte  is  exceptional  for  two  reasons.  Firstly,  there  is  no 
ionic  reactant  from  the  electrolyte.  Unlike  redox  systems,  the  system  does 
not  have  equimolar  counter -current  mass  transport  of  reactant  and  product 
ions.  Mass  transport  limitations  occur  for  dissolved  hydrogen  diffusing 
through  the  film  thickness  W,  and  for  Ht^H20  diffusing  horizontally  away 
from  the  area  of  reaction.  (However,  H2O  can  be  considered  as  a  reactant.) 
Secondly,  hydrogen  ion  has  a  high  mobility  compared  to  other  ions.  For  these 
reasons,  we  have  to  consider  anew; the  problem  of  mass  transport  polarization 
versus  ohmic  potential  gradient. 

Let  us  consider  the  mass  transport  of  product  hydrogen  ion,  with  concen¬ 
tration  denoted  by  F.  Consider  a  differential  element  dx  through  the  electrode. 
At  steady  state,  the  mass  transport  of  P  along  the  film  (assumed  linear)  is 
related  to  current  by 

1  -  -nFDp(dP/dx)(AW/q)  (52) 

A  is  the  macropore  perimeter  (per  sq.  cm.  of  electrode  face)  and  W  is  the 
mean  thickness  of  the  film,  therefore  the  cross-sectional  area  for  horizontal 
mass  transport  is  AW:q  is  the  .tortuosity  factor.  This  neglects  geometric 
effects,  but  the  geometry  of  the  pores  of  an  actual  electrode  is  likely  to  be 
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complex,  and  the  linear  approximation  is  probably  as  valid  as  assuming  some 
other  geometry,  e.g.,  a  cylindrical  geometry.  The  ohmic  effect,  again  taking 
a  linear  approximation,  is 

i  =  (l/p')(dT)/dx)(AW/q)  (53) 

Combining  and  integrating  equations  32  and  S3 

Tig-Ti  -  p'nFDp(P-Pg)  (54) 

where  T)  ^  P.  refer  to  the  polarization  and  concentration  at  the  junction  of 
'S  s 

the  thin  film  and  the  bulk  pore  electrolyte.  As  the  current  density  becomes 

greater,  hydrogen  ion  builds  up  in  the  film  and  the  specific  resistance  goes 

down.  The  minimum  specific  resistance  is  about  2  ohm  cm  and  p'  will  normally 

-4 

be  higher  than  this;  D  for  hydrogen  ion  is  about  10  ,  therefore, 

Tlg-Tl  >  20(P=^^) 

P  can  easily  increase  from  1  N  to  10  N,  giving  an  ohmic  voltage  change  of  at 
least  200  millivolts.  Thus,  in  this  case,  the  large  effective  diffusion 
coefficient  of  hydrogen  ion  combined  with  the  large  possible  concentrations 
make  ohmic  effects  appreciable  compared  to  mass  transport  effects. 

Since  there  is  no  counter -current  flow  of  reactant  product  ions,  the 
outward  mass  transport  of  hydrogen  ion  leads  to,  at  steady  state,  a  diffusive 
inflow  of  water.  Since  water  is  not  discharged,  this  must  be  matched  by  an 
outward  bulk  flow  of  electrolyte.  As  before  (see  section  3.2),  when  P  is  not 
very  high  the  outward  flow  of  water  will  not  much  increase  the  effective 
diffusion  coefficient  of  F.  We  will  assume  as  a  first  approximation  that  the 
activity  of  water  remains  1.  We  will  also  assume  that  the  specific  resistance 
of  the  film  is  constant. 

Case  1.  Reversible  Electrochemical  Reaction. 

For  a  good  surface  catalyst  such  as  highly  active  platinum  it  is 
possible  that  the  electrochemical  reaction  remains  close  to  reversible  at 
all  current  densities.  Then  the  polarization  at  any  point  x  in  the  thin 
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fllm  Is  related  Co  the  concentration  of  dissolved  hydrogen  at  the  active 
surface^  H  say^  and  the  concentration  of  such  that 

T)  -  (RT/F)  ln[(H^/H)^^^(P/Pj^)] 

Hg  Is  dissolved  H2  concentration  corresponding  to  saturation  at  the  pressure 
of  the  gas  In  the  pore.  Or 

Let  us  assume,  again  as  a  first  approximation,  that  for  concentrated  acid 
solutions  the  variation  of  P/P^  along  the  film  can  be  neglected.  Also  we 
will  assume  that  mass  transport  of  dissolved  hydrogen  occurs  across  W,  that 
Is,  It  Is  a  linear  process  perpendicular  to  the  active  surface.  It  Is  also 
assumed  that  mass  transport  resistance  In  the  gas  phase  . Is  negligible.  Thep, 
since  Adx  Is  cross  section  for  vertical  mass  transport. 


dl  -  (2FD/W)(Hj^-H)Adx 
or  H/Hjj  -  1  -  (W/2FDAHg)dl/dx 
Substituting  Into  the  previous  equation  for  polarization. 


■2Tl/b  ^ 


(W/2FnAHg)di/dx 


But  from  equation  53 


di/dx  =  (AW/qp')  d^Tl/dx^, 


therefore. 


d^Ti/dx^  -  k(l-e"^^*’) 


where 


1/k  -  W  /2FDH,-qp' 

D 


The  total  current  density  Is  given  by 


1  -  (AW/p'q)(dll/dx) 
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t  is  the  length  of  the  thin  film  up  to  the  bulk  pore  electrolyte. 

The  partial  solution  of  equation  56  can  be  obtained  as  follows.  Let 
dT|/dx  -  y,  then 

d^Ti/dx^  -  g  -  £  -  ydy/dTi 

Thus  equation  56  goes  to 

ydy/d'q  ■  k(l-e 
Separating  and  integrating^ 

y^/2  =  kT)  +  (bk/2)e‘^^^'’  +  C 

dT)/dx  =  (2kTi  +  (4bk)e"^^'^‘’  + 

From  the  boundary  condition  that  dT)/dx  “  0  at  x  «  0  (since  i  “  0  at  x  =  0] 
C  ■  -(2k1l^  +  4bke  ^o^^)  and 

dT)/dx  =  2k  [(TI'T)^)  +  (b/2)(e  -  e 


This,  of  course,  is  only  a  partial  solution  since  is  unknown. 

Let  us  consider  low  polarization  conditions  where  T)  is  low.  Then, 

2 

using  the  exponential  expansion  and  neglecting  higher  power  terms  than  T)  , 

dTi/dx  =  N/2k  (T)-Ti^  +  2  '  b  ■  2  ■  b^ 

“  (Tl^-%^)(2k/b)^/^ 

Separating  and  integrating. 


ln[(T]  /C]  =  ^  2k /b  x 


When  x  ■  0,  T]  ■  T)  ,  therefore  C  =  T)  and 


1 


ri  e 
*0 


r2  V  2  -  V  2  2*^2k/b  x  .  2  .  ^  „'/2k/b  x 


(59) 


and 


X 
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or 

T,  =  (l/2)n„  <e  *  +  e-  ’■)  '  ■ 

At  X  ■  t,  T]  T^s,  therefore 

-  +  e-''^  ’‘)/(e‘^  '  + 

*  (60) 

Therefore,  from  equation  58 

i  -  (AW/|5'q)Tig  N/2k/b  tanh(6/2k7b  t)  (61) 

The  result  obtained  is  that  i  is  directly  proportional  to  T)  ,  which  is  the 

s 

measured  polarization  (allowing,  of  course,  for  ohmic  drop  in  the  bulk  of 
the  electrolyte) .  The  slope  is 

(AW/p!,)  :  , 

-  n/(2AFDHjj)  (A/p'q)(2/b)  tanh  -y  (2FDH^^)  (2qp7b)  (t/W)  ^  (61a) 

Good  estimates  can  be  made  of  all  of  these  terms  except  t/V.  t/\f  can  be 
replaced  by  a  measurable  quantity  by  expressing  it  in  terms  of  the  limiting 
current  density  i^^. 

The  limiting  current  is  obtained  when  H  0  along  all  of  the  thin  film, 
therefore  from  equation  55, 

(62) 

/• 

atmosphere  of  H2  is  about 


=  (tA^)  milliamps/sq.  cm. 

Thus,  values  of  t/V  of  the  order  of  100  appear  reasonable,  assuming  that 
the  limiting  currents  observed  experimentally  are  due  to  the  physical  model 
discussed  above. 


ij^  =  (tA/)2FDAH^j 


For  example,  2FDH.  at  room  temperature  and  1 

“6  ^  323 

10  .  For  a  value  of  A  of  10^  cm  /cm 


-3 

i.  =  (t/W)10  amps/sq.  cm. 

1j 
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When  the  term  *7  ( 8FDHj^)  ( qp '  /b)  ( t /V)  ^  is  small  (<  0.1)  the  tanh  term 
reduces  in  the  form  tanh  x  =  x.  Then 


i  =  2(2AFDH, /b)(t/W)  (61a) 

D  *8 

or  i/i  =  2r\  /b 

L*  *S 

The  slope  is  now  independent  of  the  specific  resistance  of  the  electrolyte. 
This  is  the  case  of  pure  concentration  polarization  with  negligible  ohmic 
effect.  When  the  tanh  term  reduces  in  the  form  tanh  x  =  1,  (x  >  2) 


i  =  n/  (2AFDH^)  (A/p  '  /q)  (2/b)  -Tlg 


(62a) 


The  apparent  conductivity  is  now  proportional  to  the  square  root  of  the  true 
specific  conductivity. 

Returning  to  equation  59  to  consider  the  full  form  of  the  equation,  it 
can  be  seen  that  it  does  not  have  a  convenient  general  solution.  However, 
physical  reasohing  tells  us  that  if  the  ohmic  effect  through  the  thin  film 
is  large,  then  T)^  will  be  small.  In  fact,  at  any  given  current  density  the 
value  of  T)^  is  less  than  which  is  less  than  where  Tj^  is  the  polari¬ 
zation  which  would  be  obtained  from  simple  concentration  polarization  (of 
the  hydrogen  through  the  thin  film)  .  T]^  is  given  by 

T)^  =  (2-3)  (b/2)  log(Y:^) 

L 

Thus  T)^  is  not  large  until  i  approaches  i^^.  Since  is  much  less  than 
when  the  ohmic  effect  is  large  (that  is,  is  large)  then  T)^  will  approach 

zero  except  for  conditions  near  the  limiting  current.  Therefore,  equation 
59  can  be  put  in  the  form 


dTi/dx 

S 

From  equation  58,  then. 
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=  -/(2AFDHg)(2A/qp')  ]l/2 

As  expected,  this  reduces  to  equation  62a  when  is  small.  The  equation 
can  be  put  into  a  form  suitable  for  computation  by  noting  that  2AFDHjjt/W  =  i^ 
and  qp't/AW  =  p,  the  resistance  of  the  electrolyte  (for  1  sq.  cm.  of  electrode) 
in  the  thin  film  of  cross-section  AW  and  length  t.  Then  equation  63  can  be  put 
as 

i  =^(iJP)(h/2)  [(2Tig/b)  -  (l-e‘^"ls/b)]l/2 

Letting  2  p/b  =  r,  T)  /b  = 

s 

i  =^i^  [(2^)  -  (l-e'^S]^/2  (63a) 

For  a  given  system  i, ,  r  are  constants  while  is  the  polarization  as  multiples 
of  b(=RT/F  =  0.025  volt  at  room  temperature) .  The  form  of  the  equation  is 
shown  in  Figure  5.  It  must  be  remembered  that  the  equation  is  not  valid  at 
the  limiting  current,  therefore,  Figure  5  does  not  show  any  limiting  current. 
For  values  of  2r\  greater  than  about  2 

i  oc^y(2Tl/b)-l  (63a) 

2 

Thus  i  versus  T)  is  a  straight  line. 

Figure  6  shows  the  very  pronounced  effect  which  r  has  for  a  given  limiting 
current.  However,  it  is  unlikely  that  a  process  will  have  a  high  limiting 
current  fdi  t/W^with  a  low  v,alue  of  r(sc;  tA?).  Refering  to  . 
equation  63,  it  will  be  seen  that  parameters  which  decrease  polarization  at 
a  given  current  are  A  and  H„.  However,  if  A  is  made  too  large,  the  free  gas 

D 

space  will  become  small  and  gas  transport  restrictions  will  occur.  H^  can  be 

D 

increased  by  increasing  pressure,  according  to  a  Henry's  law  relation. 

Case  2.  Low  Exchange  Current  Reaction. 

In  the  case  considered  previously,  the  ohmic  effect  along  the  thin 
film  of  electrolyte  came  into  balance  with  the  mass  transport  of  hydrogen 
through  the  thin  film.  For  low  exchange  current  reactions,  the  ohmic  effect 
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can  also  come  into  balance  with  the  kinetic  polarization.  For  long  thin  films, 
it  is  possible  that  mass  transport  through  the  film  is  rapid  even  at  fairly 
high  currents  and  ohmic  voltage  gradients  in  the  film.  Then,  providing  the 
current  is  not  of  the  same  order  of  magnitude  as  the  limiting  current,  we 
may  set  H/H^  constant.  The  situation  now  corresponds  to  that  discussed  in 
Case  3,  section  3.4.2.  Case  3B,  that  of  high  exchange  currents,  will  clearly 
tend  to  the  reversible  condition  treated  above  as  Case  1,  section  4.2.  There¬ 
fore,  we  would  not  expect  it  to  be  possible  to  consider  it  without  allowing 
for  mass  transport  of  hydrogen  through  the  thin  film.  However,  for  Case  3A, 
that  of  low  exchange  current,  we  may  be  able  to  neglect  mass  transport  effects 
providing  we  are  well  away  from  the  limiting  current.  Then 

dl/dx  .  l^S 

and 

di/dx  =  (AW/qp')d^T^/dx^ 

Therefore, 

d^T)/dx^  =  (qp'i^S/AW)  (64) 


The  assumption  implicit  in  equation  64  is  that  the  exchange  current  is  so 
low  that  the  polarization  at  x  =  0,  T]  say,  must  be  large  (large  enough  for 
e  'o  to  be  negligible)  before  enough  current  is  produced  to  give  an 

appreciable  ohmic  effect.  As  before,  we  can  derive  an  equivalent  equation  to 
equation  43, 


d(ATi)/dx  =  [(2kb/a)  -  (2kb/a)]^'^^ 


where  k  is  now 

k  =  (qp'i^S/AW)e°'^o'''’  (64a) 

Let  a  y,  dy/dAT]  ■  (a/b)e^^^^  «  (a/b)y.  Then,  separating  and  inte¬ 

grating. 
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•J  2kb/a  X  ■ 


o  (a/b)y  "J y-1 

■  (b/a)2  tan  n/ y-1)  +  C 

-  (2b/a)  tan'^  +  C 

When  X  ■  0,  AT]  ■  0,  therefore  C  ■  0  and 

e°^Tl/b_j^  ■  (tan  x  si ka/2b)^ 
or 

(a/b)  (tIs’^q)  “  ln[(tan  t  >Jva/2b)^  +  l]  (65) 
Differentiating  with  respect  to  x, 

(a/b)dT|/dx  ■  2  ^yka/2b  tan(x  '>/ka/2b) 

Since  i  ■  (AW/qp')  (dT)/dx)^_^ 

i  =  (AW /qp ' )  n/  2kb /a  tan(t  ^  ka/2b) 

=  ^y(AWi^S/qp')(2b/cc)  tan(7  qp'i^S /AW)  (a/2b)t^  (66) 

This  is  the  explicit  relation  between  i  and  while  equation  65  allows  to 
be  calculated.  They  can  be  expressed  in  convenient  parameters  by  noting  that 
the  total  exchange  current  density  of  the  electrode  is  i^  =  while  the 

resistance  of  the  thin  film  of  electrolyte  (for  1  sq.  cm.  of  electrode  face) 
is  qp't/AW  ■  p  say,  then 

i  -  /(I  /p)(2b/a)  tanf/"(i  5)(a/2b)  (66a) 

o  o 
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For  low  values  of  the  tan,  that  is,  low  values  of  (and  current),  tan 

X  ■  X,  therefore 

i  -  i  e°^o/b 
o 

and  ‘Hg  ”  'Hq-  This  is  the  required  logical  result,  since  we  assumed  an 
irreversible  Tafel  form  as  our  basic  equation.  As  the  tangent  becomes  large, 
from  equation  65,  neglecting  the  1,  tan  =  exi((a/2b(Ti  -l)  )]  and 

8  O 

i  -  f  (i^/p)(2b/a)  e°^s/2b 

This  is  the  corresponding  form  to  equation  45.  The  full  form  of  the  equation 
can  be  computed  considering  as  parameters  pQ!/2b  =  r  say,  OTj/b  ■  T],  then 

i  =  n/I  /r  tan(  JTT  e'^o^^) 

o  o 

Tlg-T^^  *  ln[ (tan  J i^r  +  l] 

or 

i/i^  =  (l/A)e^o^^  tan(Ae^o'^^) 

^s"\  “  ^  ^  log[^an(Ae^o^^))^  +  l] 
where  A  =  JTr  (dimensionless) 

When  's/  Igr  e^°^^  =  ir/2,  the  tangent  is  infinity  and  T)^  is  infinity.  There¬ 
fore  the  maximum  value  of  for  a  given  system  is  given  by 

(Ti  )  =2  ln(7r/2»yTT)  (67) 

'o  max  o 

(ti  )  cannot  be  small,  since  we  have  assumed  an  irreversible  reaction, 

•o  max  ’  ’ 

therefore  equation  67  only  applies  when  T)^  is  greater  than  about  2. 

The  forms  of  the  results  are  shown  in  Figure  7,  where  i/i^  is  plotted 
versus  T]^,  T)^  for  various  values  of  i^r.  [The  order  of  magnitude  of  r  can  be 
decided  by  taking  qp'  as  of  the  order  of  10,  A  as  of  the  order  of  10^,  tA7  as 
in  the  range  10  to  1000  and  a/2b  as  approximately  10.  Then  r  would  be  expected 
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to  range  from  1  to  100.]  As  expected,  for  a  given  value  of  a  decrease 
in  the  value  of  r  gives  a  lower  polarization  at  high  current  densities  and 
a  more  extended  true  Tafel  line.  The  current  density  at  which  the  double 
Tafel  slope  line  takes  over  is  higher.  Alternatively,  for  a  given  value  of 
r  a  bigger  exchange  current  has  the  same  effect. 

4.3.  Thin  Film  Model;  Hydrogen  Fuel  in  Alkaline  Electrolyte 
In  strongly  alkaline  electrolyte  the  reaction  probably 
involves  hydroxyl  (14)^.  as  follows, 

[h]  +  OH'  ^  H2O  +  e 

In  this  case  we  do  have  an  ionic  reactant  from  the  electrolyte  and 
therefore  we  would  expect  mass  transport  of  OH'  in  the  thin  film  to  be  a 
critical  factor.  As  in  the  case  of  redox  reactions,  the  concentration  of 
the  hydroxyl  ion  will  tend  to  zero  at  x  =  0  and  equation  33  (and  its  reduced 
forms)  may  apply.  It  should  be  noted  that  the  thin  film  of  electrolyte  may 
have  a  low  conductivity  due  to  its  narrow  cross-section.  However,  the 
effective  mass  transport  of  ion  in  the  thin  film  will  also  be  low,  by  exactly 
the  same  ratio.  Therefore,  what  was  said  before  about  the  relative  effects 
of  mass  transport  to  ohmic  effects  will  still  apply  (see  section  3.4.1). 

The  ohmic  effect  should  be  considered  as  a  disturbance  of  the  basic  relation 
of  equation  33,  and  it  will  give  somewhat  higher  total  polarizations  than 
predicted  by  equation  33.  The  effect  will  be  most  marked  at  high  current 
densities,  therefore  it  leads  to  a  more  gradual  approach  to  the  limiting 
current  than  expected. 

The  relative  ease  of  mass  transport  at  the  exterior  to  that  in  the 
Interior,  Is  given  for  the  thin  film  model,  by 

if/  =  qt/8AW  (68) 

Comparing  with  equation  32,  it  is  seen  that  e  is  replaced  with  AW.  AW 

3 

cannot  be  greater  than  €.  For  A  =  10  ,  W  =  1  micron,  AW  is  1/10.  This 
compares  with  a  value  of  €  of  about  1/3.  Of  course,  W  might  be  considerably 
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less  than  1  micron. 

4.4.  Non  Thln-fllm  Model;  Saturated  Electrode 

Another  physical  model  which  might  be  applicable  under 
certain  conditions  is  one  where  no  extended  thin  film  exists.  This  would 
apply  to  a  fully  saturated  electrode.  In  this  case  the  reactant  gas  has 
to  diffuse  into  the  electrol3'te  in  the  macropores.  When  current  is  drawn, 
the  reactant  gas  in  the  interior  will  be  rapidly  used  up  and  the  reaction 
will  be  concentrated  at  the  gas -electrolyte  face.  Ohmic  polarization 
increases  from  this  face  towards  the  bulk  electrolyte  and  this  will  have 
the  effect  of  speeding  up  internal  reaction,  giving  rise  to  low  internal 
dissolved  gas  concentrations,  which  again  means  that  the  reaction  is  concen¬ 
trated  towards  the  gas.  entry  face.  This  is  opposite  to  the  ionic  mass 
transport  case,  in  which  the  concentration  of  reactant  ion  decreases  from 
the  free  electrolyte  face  back  into  the  electrode.  For  the  ionic  case, 
reaction  is  concentrated  at  the  electrode-bulk  electrolyte  face  and  conse¬ 
quently  the  path  of  the  current-carrying  ions  does  not  penetrate  deeply 
into  the  electrode.  For  the  dissolved  gas  case,  the  current-carrying  ions 
have  to  travel  right  through  the  pores  of  the  electrode  to  (or  from)  the 
reaction  site  and,  therefore,  ohmic  effects  within  the  pore  electrolyte  are 
much  greater  than  for  the  ionic  case. 

Let  us  consider  the  case  where  the  electrochemical  reaction  is  reversible, 
where  the  principal  mass  transport  effect  is  that  of  dissolved  reactant  gas 
and  where  the  dissolved  gas  has  to  diffuse  linearly  along  the  axis  of  the 
pore.  If  the  concentration  of  dissolved  gas  is  R,  and  the  equilibrium 
concentration  at  the  gas  pressure  used  is  R^, 

1)  =  (RT/nF)ln(R^^/R) 

If  D  is  the  effective  diffusion  coefficient  of  dissolved  gas, 

di/dx  -  nFD(d\/dx^)  (69) 


Also 


i  ■  (1/p)  dTl/dx 
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The  boundary  conditions  are  i  =  0  at  x  =  0,  R  =  at  x  =  0  and  i^  * 
-nFD(dR/dx)^_Q,  where  i^  is  the  final  current  density.  Integrating 
equation  69  from  x  =  0  to  x  =  x, 

i  =  nFD(dR/dx)  -  nFD(dR/dx)^Q  =  i^.  +  nFD  dR/dx 

Therefore, 

dTl  =  pi^x  +  pnFD  dR 

integrating  from  x  =  0,  ri  *  to  x  =  t, 

Tl-T)^  «  pnFD(R-R^)  +  ij.px 
or 

R/R^  =  1  -  (ij.x/nFDR^)  +  (T)-T)^) /(pnFDR^) 

The  measured  polarization,  at  the  electrolyte  face  of  the  electrode, 
X  ■  t,  therefore 

e-”’!'’'  .  (R^R^XR/R^) 

-  (RQ/Rb)[l  -  (i,./nFDR^)  +  (Tl-T)^)/(pnFDR^)] 
or 

i  =  (nFDR^^/t)[(R^/R,^)-e'"^^^]  +  (Tl-Tl^) /tp 

If  there  is  no  gas  mass  transport  limitation  to  x  =  0  then  R^/Rj^  “  1 

Tl  ■  0.  Then 
'o 

i  -  (nFDR^^/t)(l  -  +  T^/tp 

When  T]  is  small, 

i  =  'n[(l/tp)  +  (n^FDR^/bt)] 

-  (l/tp)Tl[l  +  (pnWjj/b)] 


(70) 
is  at 


(71) 

and 

(71a) 


(71b) 
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Thus  T]  is  linearly  dependent  on  i.  tp  is  the  ohmic  resistance  of  the 

electrolyte  in  the  pores,  therefore  the  effective  resistance  is  less  as 

the  mass  transport  of  dissolved  gas  is  greater.  For  strongly  conducting 

-2 

electrolytes  the  factor  pnFDRj^/b  is  of  the  order  of  10  M,  where  M  is  the 
molar  concentration  of  dissolved  gas.  Since  M  is  low,  the  term  containing 
is  normally  negligible.  The  physical  meaning  of  this  is  that  the  dis¬ 
solved  gas  mass  transport  is  so  slow  that  the  reaction  occurs  in  a  thin 
layer  near  x  =  0.  Then  all  the  current  has  to  traverse  the  complete  - 
thickness  of  the  electrode  and  the  ohmic  voltage  drop  is  the  maximum 
possible. 

When  T)  is  large,  equation  71  goes  to 

i  =  (nFDR^/t)  +  (ll-Tl^)/tp 

For  dissolved  gases  and  strongly  conducting  electrolytes  the  term  in  is 
negligible  and 

i(tp)  =  TJ-Ti^  (72) 

This,  of  course,  is  the  expected  Ohm's  law  relation.  The  total  polarization 
is  made  up  of  an  ohmic  part,  T)-T)^,  and  a  mass  transport  part  T)^.  The  mass 
transport  part  may  consist  of  diffusion  through  a  thin  film  of  electrolyte  to 
reach  reaction  surface  at  x  =  0.  If  so,  the  limiting  current  will  be  small. 
If  no  thin  film  exists,  then  electrolyte  at  x  =  0  is  in  equilibrium  with  the 
reactant  gas  pressure  at  the  electrolyte  surface  and  the  mass  transport 
limitation  will  be  that  of  gas  transport  to  the  surface.  However,  since  the 
reaction  is  concentrated  more  and  more  into  a  thin  layer  at  x  -  0  as  ohmic 
polarization  increases,  the  true  current  density  becomes  high  and  activation 
polarization  may  be  encountered. 

4.5.  Surface  Effects. 

When  a  catalyst  is  impregnated  or  plated  onto  the  internal 
surface  of  a  porous  electrode,  its  principal  effect  is  to  Increase  the 
exchange  current  of  the  reaction  involved,  which  forces  the  system  towards 


-43- 


a  pure  mass  transport -ohmic  polarization  condition.  This  condition  represents 
the  best  possible  performance  obtainable  with  the  given  electrode  structure. 
Once  this  condition  is  reached,  improved  catalysis  of  the  reaction  will  have 
no  further  benefit.  Electrode  performance  can  then  be  improved  by  making 
the  electrode  thinner  and  more  open,  thereby  reducing  internal  mass  transport 
and  ohmic  limitations.  But,  by  doing  this  the  effective  catalytic  area  is 
reduced  so  that  activation  polarization  may  again  come  into  play.  For  the 
thin  film  model  of  a  gas  diffusion  electrode,  reducing  the  thickness  of  the 
electrode  beyond  a  certain  limit  may  reduce  t,  the  length  of  the  thin  film. 
Polarization  at  low  currents  is  then  reduced,  but  the  limiting  current 
becomes  less.  Clearly,  design  of  the  optimum  electrode  requires  the  balancing 
of  a  number  of  different  factors.  At  present,  the  theories  described  above 
have  not  been  sufficiently  investigated  experimentally  to  show  if  they  are 
correct.  Also,  the  important  parameters,  such  as  t/W,  have  not  been  deter¬ 
mined. 

When  a  catalytic  surface  is  micro -roughened,  the  principal  effect  is 
again  that  of  increasing  the  effective  exchange  current.  A  secondary  effect 
may  be  that  the  roughened  surface  holds  electrolyte,  thus  giving  a  greater 
cross  section  of  electrolyte  for  ionic  mass  transport  and  ionic  conduction. 

The  gas  diffusion  limiting  current  would  still  be  determined  by  t/yi,  but 
ionic  mass  transport  and  ohmic  polarization  effects  would  be  smaller.  If 
this  were  true,  it  might  be  possible  for  W  to  be  very  thin,  giving  large 
limiting  currents,  without  the  ohmic  resistance,  for  example,  tending  to 
large  values.  This  would  appear  to  be  a  more  refined  model  well  worth 
investigating . 

Several  surface  chemistry  effects  can  disturb  some  of  the  analysis 
described.  Firstly,  some  of  the  work  is  based  on  the  applicability  of  a 
simple  Tafel  equation  for  the  reaction  studied.  If  several  steps  in  the 
over-all  reaction  proceed  at  comparable  rate,  then  the  rate  equation  is 
more  complex  and  the  form  of  the  rate  equation  may  change  with  current 
density.  However,  it  is  often  found  that  the  Tafel  equation  applies  over 
wide  ranges  of  current  density.  Secondly,  the  Tafel  equation  does  not 
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always  take  into  account  the  variation  of  the  amount  of  surface  complex  with 
conditions.  Since  the  electrode  reactions  are  surface  reactions,  the 
activities  of  covered  surface  and  bare  surface  should  appear  in  the  rate 
equations.  Where  gas  adsorption-desorption  is  involved,  the  change  in 
surface  coverage  can  lead  to  limiting  currents  of  chemisorption  or  de¬ 
sorption.  Since  the  Tafel  equation  is  found  so  often,  it  might  be  postu¬ 
lated  that  the  complex,  for  redox  and  gas  reactions,  has  a  Temkin  isotherm 
for  its  surface  coverage  versus  heat  of  reaction  relation.  In  this  case, 
large  changes  in  current  density  can  occur  with  only  small  changes  in 
surface  coverage  (14).. 

(5)  UMSTEADY  STATE  REDOX  SYSTEMS..  • 

5.1.  Introduction. 

The  unsteady  state  relations  for  porous  electrodes  have  not 
received  much  study.  It  would  appear  reasonable  to  obtain  experimental 
confirmation  of  the  simpler  steady  state  relations  before  spending  much 
time  on  the  unsteady  state  theory.  For  a  redox  system,  as  in  Figure  1, 
we  can  assume  as  a  first  approximation  that  ohmic  gradients  within  the 
pore  electrolyte  are  small. 

Considering  an  element  dx,  the  rate  of  build-up  of  reactant  is 
€(SiR/Sit)dx  (t  being  time).  The  sum  rate  of  diffusion  of  reactant  into 
the  element  is  D(S^/9x^)dx,  where  D  is  understood  to  be  an  effective  mass 
transport  factor.  The  rate  of  conversion  of  R  to  P  in  the  element  is  given 
by 

di^/nF  =  (i^S/nF)[(R/Rj^)e“^^’^-(P/Pj^)e'^^'^^’^^'’]dx 
ij^  is  the  faradalc  furrent.  Then 

eSR/6t  =  D  x^-(i^S/nF)[(R/R^^)e°^/'’-(P/Pj^)e'^^'“)^/*’]  (73) 

The  current  density  at  time  t  is  given  by 
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i  =  i-  +  1 
f  c 

d 

=  /  di^  +  CdTl/dt 
o 

d  is  the  thickness  of  the  electrode,  i^  is  the  capacitative  current,  C  is 
the  capacity.  The  boundary  conditions  are  (Sr/Sx)  =  0  at  all  values 
of  t,  and  D(3r/Sx)^_^  =  D^(5R/5y)y_^,  where  is  the  effective  mass 
transport  factor  in  the  free  electrolyte  and  y  is  distance  into  the 
electrolyte  from  the  electrode-electrolyte  Interface. 

Clearly,  the  complete  solution  of  these  equations  for  a  number  of 
initial  conditions  is  not  easy.  However,  two  simplified  cases  can  be 
readily  examined. 

5.2.  Large  Applied  Current 

Consider  a  redox  electrode  at  open  circuit  with  uniform 
concentrations  throughout.  Let  a  large  constant  current  (greater  than  the 
limiting  current)  be  applied  to  the  electrode.  For  short  periods  of  time, 
the  fresh  reactant  diffused  in  from  the  exterior  will  be  small  and  the 
current  will  be  supported  mainly  by  the  existing  reactant  in  the  interior 
of  the  electrode.  The  system  is  as  illustrated  in  Figure  8.  Let  us 
assume  that  the  concentration  of  reactant  in  the  interior  is  constant 
with  X  and  thus  goes  down  as  a  horizontal  line  in  Figure  8.  Assuming 
negligible  reactant  transfer  into  the  electrode  from  the  outside. 


i  =  CdT|/dt  +  i^  =■  constant  applied  current 
€  dR/dt  =  -(i^S/nF)[(R/R^^)e®'^'^^-(P/P^^)e’^^‘°‘^''^^*’] 
ij  -  i^Sd[(R/R^^)e“^^'’-(P/P^^)e'*'^'“^’l^*’] 

Then, 

dR/dt  =  -(i  -  CdT)/dt) /enFd  (74) 


Integrating  from  t  ■  0,  R 


Rb^  T|  -  0  to  t  -  t,  R  -  R,  1)  ■  n. 
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Rj^-R  =  (it/denF)-(CTl/denF) 

(75) 

R/R^  -  1  -  it/denFR^^  +  CTl/denFRj^ 

(75a) 

Folarlzatlon  becomes  large  when  «  1  and  a  "transition  time"  t 

can  be  defined  by  this  condition.  Then 


T  =  (d€nFRj^/i)  +  (CT)/i) 

(76) 

Then  R/R^^  =  (i/d  nFRj^)  (T-t) 

(75b) 

For  example^  consider  d  =  0.1  cm.,  e  =  1/3,  n  =  1,  Rj^  ■  0.1  M, 

T  2i  0.3/i  +  Cn/i 

For  a  charge  current  density  of  0.1  amps  per  sq.  cm.  and  a  capacity  of 
4  5 

10  pF  per  sq.  cm.  (=  10  pF  per  cubic  cm.  in  this  case),  and  assuming  that 
the  transition  time  is  reached  by  the  time  T)  ■  1/2  volt, 

T  ^  3  +  1/20  secs. 

Under  these  conditions  the  effect  of  double  layer  capacity  on  the  transition 
time  is  negligible.  In  3  seconds,  mass  transport  from  outside  the  electrode 
will  be  fairly  small. 

The  voltage-time  relation  is  obtained  by  substituting  for  R/R^  from 
equation  75a  to  give 

i  =  CdTj/dt  +  i^Sd[(i/d€nFRj^)(T-t)e®^^^-(P/P^j)e'^^"“^'^^'’]  (77) 

This  cannot  be  usefully  integrated  over  the  entire  range.  However,  some 
limiting  cases  can  be  considered. 

Case  1.  Low  Exchantie  Current  Reaction,  Irreversible  Conditions. 
Consider  the  period  of  time  when  the  capacitative  current  is  small 
and  Cdr|/dt  is  small.  This  applies  when  the  faradaic  current  is  small  until 
a  fair  amount  of  polarization  is  present;  the  initial  current  is  then  mainly 
capacitative  until  this  polarization  is  reached.  When  the  reverse  reaction 
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term  is  negligible^ 

i  2;  iQSd[(i/denFR^)(T-t)]e‘^^’^ 

Replacing  d€nFR|^/i  by  T(when  CT)/i  is  negligible  near  the  transition  time) 

-r\  2;  (b/a)ln(i^Sd/i)  +  (b/a)ln(l-t/T)  (78) 

Thus  a  plot  of  T)  versus  log(l-t/T)  is  a  straight  line  and  an  can  be  obtained 
from  the  slope  and  i^Sd  from  the  intercept,  [if  C  is  known,  t  can  be  corrected 
for  CT)  and  i  for  CdT)/dt.] 

As  t  approaches  the  transition  time  and  the  polarization  becomes  large, 
reaction  on  the  external  surface  will  start  to  carry  current.  This  in  turn 
will  have  a  transition  time,  defined  in  the  usual  way  by 

T  =  (Rjj/2p^(nF)  V 

2 

Under  the  conditions  used  before,  that  is,  R,  =0.1M,  i=0.1  amps /cm.  , 

-52  “ 

and  assuming  D  to  be  about  10  cm.  /sec.,  t  ^  0.075  seconds.  Thus  the 

extension  of  the  total  transition  time  due  to  the  external  current  is  small 

under  these  circumstances.  However,  the  second  transition  time  is  proportion- 
2 

al  to  (Rj^/i)  ,  whereas  the  first  is  proportional  to  R^^/i.  If  made 

too  large,  then  the  disturbance  of  the  transition  time  of  equation  76  by  the 
external  face  current  may  be  excessive.  It  can  be  made  negligible  if  d,  the 
thickness  of  the  electrode,  can  be  varied  to  suit  our  purposes. 

An  estimate  of  the  effect  of  mass  transport  into  the  interior  can  be 
made  as  follows.  Consider,  as  before,  that  the  concentration  R  in  the  porous 
electrode  is  constant  up  to  the  face.  Then  the  rate  of  mass  transport  into 
the  Interior  is,  in  the  usual  manner  ^  (19),  given  by 

dQ/dt  =  (R. -R)(D/Tr)^^^/t^^^  Km.  moles  (79) 

^  2 - 

cm.  sec. 

Q  is  the  total  amount  transported  in  time  t.  The  total  mass  transport  in 
time  T  is 

/This  corresponds  to  the  linear,  serai-infinite  heat  transfer  case  of  a  block 
of  material  initially  at  uniform  temperature  R^,  with  a  fixed  temperature 
of  R  applied  to  one  face. 
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Q  =  2(R^^-R)(D/Tr)^^^  (80) 

But  R  is  changing  with  time  as  given  approximately  by  equation  75  and 
applying  the  principle  of  superposition, 

Q  =  /  2(DyTr)^^^  (T-t)^^^  d(R  -R)  (81) 

o  “ 

Substituting  for  Rjj"R  from  equation  75  (neglecting  the  capacity  term) 

Q  =  /2;(p/Tr)^^^(T-t)’-^^(i/denF)dt 
o 

=  (4/3).(D/7r)^^^(i/denF)T^^^  (82) 

The  total  charge  corresponding  to  Q  is  nFQ,  while  the  charge  corresponding 
to  i  is  iT.  Therefore, 

nFQ/iT  =  (4/3)(D/Tr)^^^(l/de)T’-^^  (83) 

Taking  D  as  10  d  as  0.1  cm.,  €  as  1/3  and  T  as  3  seconds, 
nFQ/ix  2;  12%. 

Thus  a  nominal  transition  time  of  3  seconds  would  actually  give  rise  to  a 
transition  time  about  10%  in  excess  of  three  seconds.  This  predicted 
disturbance  is  probably  too  high,  since  in  practice  part  of  the  concentration 
gradient  occurs  within  the  solid  rather  than  within  the  free  electrolyte,  and 
the  effective  diffusion  coefficient  within  the  solid  is  1/5  to  1/10  of  the 
free  value. 

Case  2.  Reversible  Reaction. 

For  a  reversible  reaction  the  voltage -concentration  relation  is 
T1  =  (b)ln[(R,^/R)(P/P^^)3 

Since  we  are  assuming  negligible  diffusion  in  or  out  of  the  electrode  R  +  P  > 

R,  +  P,  and 
D  0 
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T)  =  (b)ln[(R^^/R)(Rj^  +  -  R)/?^] 

or,:  if  Rjj/P^  -  7, 

-T1  =  (b)ln[(R/Rj^)/(l  +  7  -  -jR/R^)] 

Substituting  from  equation  75b,  assuming  the  capacity  term  in  equation  76 
is  small, 

-Ti  -  (b)ln[(l-t/T) /(I  +  7-7,  (l-t/x))'}  (84) 

For  7  not  large  and  t  -»  t, 

-T)  2;  (b)ln[(l-t/T)/(l  +7)]  (84a) 

Again,  a  plot  of  log(l-t/T)  versus  should  be  a  straight  line.  In  this 
case  t  is  zero  where  is  nearly  zero,  whereas  in  equation  78,  the  extra¬ 
polation  to  t  =  0  does  not  give  a  value  of  T)  near  zero. 

5.3.  Current  Interruption. 

Consider  a  redox  electrode  operating  at  a  steady  current, 
with  negligible  internal  ohmic  polarization.  When  the  current  is  shut  off 
suddenly,  the  electrochemical  reaction  will  continue  and  will  charge  (or 
discharge)  the  electrical  double  layer  until  reversible  conditions  are 
obtained.  The  relations  are  now 

i  +  i,,  =  0 
c  f 

i^  =  C(dT)/dt) 

Case  1.  Irreversible  Reaction. 

The  relation  of  faradaic  current  to  polarization  is,  as  before, 
di^  =  i^S[(R/R^)e°^^’^]dx 

For  negligible  internal  ohmic  effects,  this  can  be  written  as 
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=  i^S  /  (R/R^)dx 

o 

For  a  short  period  of  time  after  current  shut-off,  the  values  of  R/R^  will 
be  essentially  constant  due  to  the  slowness  of  mass  transport.  Then, 

CSn/St  +  i^S  ° 

Separating  and  integrating,  with  for  the  steady  current  polarization, 

s 

(i^S/(R/Rjjd3(a/Cb)t  .  ea(Tl8-Tl)/b  _  ^ 


When  is  large. 


-■q  “  (b/a)lnt  +  (b/a)ln[(a/Cb)i^S/(R/Rj^)dx]  (86) 


A  plot  of  q  versus  log  t  is  a  straight  line.  If  no  concentration  gradients 
exist  then  equation  86  is  the  usual  equation  (14),  with  i^^Sd  as  the  exchange 
current.  If  R/Rj^  does  vary,  then  it  is  convenient  to  introduce  the  pre- 
interruption  current  i, 

i  -  1  S  e®^s'^’^/(R/R  )dx 

O  Q  D 

Then 


q^-q  =  (b/a)lnt  +  (b/a)ln[(a/Cb)i]  (86a) 


The  extrapolation  to  q^-q  =  0  gives 

t  =  Cb/ai  (87) 

o 


C  can  be  estimated  if  nd!  is  known. 

The  equations  can  only  apply  at  short  times,  before  appreciable  mass 
transfer  occurs,  for  highly  irreversible  reactions.  It  should  be  noted 
that  any  ohmic  polarization  in  the  electrolyte  decays  with  time  according 
to  the  decay  of  i^.  If  AT)  is  ohmic  polarization,  it  is  readily  shown  that 


ln(Z»i/Aqg)  -  (a/b)(q-q^) 


(88) 


» 
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Figure  9  illustrates  the  general  shape  of  the  curve. 

Case  2.  Reversible  Reaction. 

As  discussed  previously  (see  section  3),  a  completely  reversible 
electrochemical  reaction  with  negligible  internal  ohmic  effect  has  a  uniform 
concentration  throughout  the  solid.  When  current  is  shut  off,  reactant 
diffuses  to  the  surface  and  Cries  to  decrease  the  concentration  polarization. 
But  this  change  in  polarization  drives  product  in  the  interior  of  the 
electrode  to  reactant.  An  electron  balance  requires  that 

electrons  taken  up  per  second  by  P  ->  R  in  the  interior  of  the 
electrode  -  electrons  produced  per  second  by  R  -»  P  at  surface 
of  electrode  =  CdT)/dt.  (89) 

The  situation  can  be  envisaged  as  in  Figure  10.  As  time  progress  R  increases 
and  P  decreases  towards  the  equilibrium  bulk  concentrations.  It  is  clear  from 
equation  89  that  we  cannot  consider  the  reaction  at  the  surface  to  be  zero 
when  the  current  is  Interrupted. 

The  mathematical  analysis  of  this  situation  has  not  yet  been  completed. 

(6)  GSAS  PHASE  TRANSPORT. 

6.1.  Introduction 

For  a  porous  electrode  fed  with  gas  at  one  face,  it  is  always 
possible  that  gas  concentration  gradients  exist  within  the  pores  of  the 
electrode.  Although  mass  transport  of  ions  or  dissolved  gases  in  the  electro¬ 
lyte  will  normally  be  slower  than  gas  mass  transport,  the  system  may  be  such 
that  long  gas  diffusion  paths  exist  compared  to  short  paths  in  the  electrolyte 
It  is  also  possible  that  concentration  gradients  exist  due  to  gas  transport  to 
the  face  of  the  electrode,  since  practical  fuel  cell  assemblies  will  use  thin 
electrode  jackets  which  will  probably  put  the  gas  flow  in  the  laminar  flow 
regime.  In  this  section  we  will  be  concerned  with  concentration  gradients 
within  the  porous  electrode,  between  the  gas  entry  face  and  the  electrolyte 
at  which  utilization  occurs.  Two  limiting  cases  occur,  gas  diffusion  at 
constant  pressure  and  gas  permeability  with  varying  pressure. 
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6.2  Gas  Diffusion. 

In  the  consideration  of  gas  diffusion  through  porous  media^ 
two  principal  types  of  diffusion  can  be  distinguished.  Firstly,  there  is 
normal  diffusion  in  which  one  gas  diffuses  through  another,  with  negligible 
influence  of  the  pore  walls  on  the  rate  of  diffusion.  This  applies  when 
the  mean  free  path  of  the  molecules  is  much  less  than  the  pore  diameter. 
Secondly,  when  the  mean  free  path  of  the  molecules  is  much  greater  than 
the  pore  diameter,  Knudsen  diffusion  occurs.  In  this  case,  the  resistance 
to  mass  transport  is  due  to  collisions  with  the  walls  and  the  presence  of 
other  gases  has  little  influence  on  the  diffusion  of  the  component  consider¬ 
ed.  Some  specific  cases  are  as  follows. 

Equimolar  Counter -Current  Normal  Diffusion  of  Reactant  and  Product.  Two 
Components . 

The  notation  and  presentation  follows  that  of  McCarty  and  Mason 
(20),  Let  be  the  molar  flux  (moles  per  sq.  cm.  per  second)  of  gas  1 
and  J2  the  molar  flux  of  gas  2  (in  the  same  direction) .  The  equations  of 
diffusion  for  one  dimensional  diffusion  are 

Jj^  =  -Dj^(dnj^/dx)  +  nj^v 

J2  “  -D2(<in2/dx)  +  n2V  (90) 

3 

nj^,  n2  are  the  concentrations  in  moles  per  cm.  and  v  is  the  velocity  of 
bulk  flow.  Since  the  mass  of  one  gas  diffusing  in  one  direction  does  not 
in  general  equal  the  mass  diffusing  in  the  other  direction,  a  definite  net 
mass  flow  in  one  direction  is  present.  The  velocity  of  flow  is  the  net 
mass  flow  divided  by  the  density  (n^m^^  +  n2m2  in  this  case,  where  m^^,  m2 
are  the  respective  molecular  weights) .  This  net  flow  must,  of  course,  obey 
the  laws  of  fluid  dynamics.  The  diffusion  coefficients  in  equation  90  are 
quite  arbitrary;  they  are  defined  by  equation  90,  and  they  may  be  concen¬ 
tration  dependent.  *  -Dj^dn^/dx  only  when  v  =  0,  therefore  the  D  values 
are  equl-mass  diffusion  coefficients.  In  other  words,  we  can  consider  a 
mass  flow  along  the  axis,  carrying  n^,  n2  with  it.  If  a  plane  is  moved 


-53- 


along  with  the  fluid,  at  velocity  v  then  the  diffusion  rates  of  components 
1  and  2  with  respect  to  this  plane  are  given  by  -Dj^dn^/dx  and  -D2dn2/dx 
respectively . 

To  procede  further  we  must  consider  the  fluid  dynamics  of  the  mass 
transfer.  For  one  dimensional  transport  the  equation  of  motion  is  (21) 

9p/Sx  +  Sr/Sy  +  pvSv/Sx  =  0  (91) 


p  is  pressure,  T  is  viscous  force,  p  is  density.  The  pressure  term  represents 
the  driving  force,  the  viscous  term  represents  the  viscous  resistance  and  the 
density-velocity  term  represents  the  bulk  momentum  term.  For  low  values  of 
p  and  V  and  negligible  viscous  force,  dp/dx  is  small.  If  we  assume  this  to 
be  true,  then  n^  +  n2  ■  constant  *  N  say.  For  steady  state  conditions,  then, 
Bnj^/dx  ■  -dn^/dx  and  equation  90  goes  to 

=  -Dj^dOj^/dx  +  nj^v 
J2  *  D2dnj^/dx  +  n2V 

For  equimolar  counter -current  mass  transfer  =  -J2  and 


0  -  -(Dj^-D2)dnj^/dx  +  v(nj^  +  n2) 


or 

V  =  (1/N)  (Dj^-D2)dnj^/dx 
Substituting  into  equation  90, 

“  -[(n2/N)D^  +  (n^/N)D2]dnj^/dx  (92) 

or 

Jl  =  (^3) 

Dj^2  conventional  equimolar  binary  diffusion  coefficient  of  kinetic 

theory  and  is  Independent  of  the  partial  concentration  of  the  components. 

We  see,  therefore,  that  in  this  case  the  mass  transport  of  reactant  can  be 
represented  by  a  Flck's  law  form,  equation  93.  The  effective  D^2  ® 
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porous  electrode,  based  on  unit  geometric  area,  is  smaller  than  the  true 
diffusion  coefficient  because  the  effective  diffusion  area  is  the  porosity. 

In  addition,  if  dx  is  linear  dimension,  the  actual  path  length  is  qx,  where 
q  is  a  tortuosity  factor.  Thus,  as  before, 

^°12^ effective  “  °12 

For  a  diffusion  path  of  t  from  the  entry  face  to  the  reaction  zone,  equation 
93  can  be  integrated  in  the  usual  way  to  give 

^o-nt  =  (1/Di2>e 
The  limiting  current  is 

♦ 

‘l  ■ 

Then 

n  /n  =  1  -  i/i 
to  L 

It  should  be  noted  that  this  is  an  approximation  which  applies  at  relatively 
low  rates  of  mass  transfer,  since  viscous  forces  have  been  neglected.  Negligible 
pressure  gradient  exists  in  the  electrode. 

Two  Component.  Non-Eoulmolar.  Normal  Diffusion 
Equation  90  can  be  expressed  as 

Jl  -  +  (n^/N)J  (95) 

J  >  +  J2  *  net  molar  transport.  With  the  same  assumption  as  before,  that 

negligible  pressure  gradient  exists  (therefore  N  *  constant)  even  though  net 
mass  transfer  is  occuring,  this  equation  can  be  integrated  if  the  flux  ratio 
J1/J2  is  known.  An  important  case  is  that  of  a  two  component  gas  mixture  in 
which  one  component  reacts  to  a  non-gaseous  product  and  the  other  component 
is  inert.  Air  reacting  at  an  oxygen  electrode  is  an  important  example.  In 
this  case,  at  x  ■  0  the  respective  concentrations  of  reactant  and  inert  are 
n^Q  and  n2Q,  say,  is  the  rate,  while  J2  is  zero.  Then 
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Jl  =  -D^2  dn^/dx  +  (n^/N)Jj^ 
or 

Integrating  from  x  =  0  to  x  =  t, 

Jl  =  ln[(l-n^/N)/(l-n^Q/N)] 

or 

Jl  =  (Di2/t)N  ln|(N-n^) /n2Q]  (96) 

where  ni  now  refers  to  the  concentration  at  the  reaction  site.  The  limiting 
current  is  nFJi  or 

=  (nFDj^2/'^>N  ln(N/n2Q)  (97) 

3 

N  is  in  gm.  moles  per  cm.  ,  and  at  1  atmosphere  and  room  temperature  it  is 
-4 

about  (1/2)10  .  Equation  97  can  be  put  into  a  form  comparable  to  a  simple 

Pick's  law  form  modified  by  a  correcting  factor, 

1l  =  (nFDi2/t)niQ  f  (97a) 

where 

f  =  (N/niQ)2-3  log(N/n2o) 

The  values  of  f  as  a  function  of  the  molar  fraction  of  inert  present  are 
shown  in  Table  1. 

Table  1 .  Correction  to  Pick's  Law  When  Inert  Gas  is  Present. 

Molar  fraction  of  inert  Correcting  factor 

in  bulk  gas  stream 

_ ^20^^^ _  ^ _ 

0 

0.0105 
0.0446 
0.204 
0.550 
1.285 
2.07 
4.55 
6.91 


1.0 

0.9 

0.8 

0.6 

0.4 

0.2 

0.1 

0.01 

0.001 
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For  air,  with  a  mole  fraction  of  nitrogen  of  4/5,  the  correction  factor  is 
about  1/5,  therefore  the  limiting  gas  diffusion  rate  is  1/5  of  that  expected 
from  a  simple  Pick's  law  form  with  1/5  partial  pressure  of  oxygen.  Oxygen 
with  less  than  afew%inert  impurities  would  give  an  f  of  about  5,  therefore, 
at  the  same  total  pressure  the  oxygen  would  give  a  limiting  diffusion  rate 
about  125  times  larger  than  air. 

As  the  percentage  of  inerts  is  reduced  the  limiting  diffusion  rate 
becomes  very  large,  tending  to  infinity.  However,  it  is  important  to  note 
that  this  is  a  spurious  prediction  for  flow  through  porous  materials.  The 
case  of  zero  impurity  is  the  same  as  mass  transport  of  one  gas,  the  reacting 
gas,  with  no  other  component  present.  This  clearly  becomes  a  permeability  or 
bulk  flow  system  when  the  rate  of  flow  becomes  high.  Thus  an  upper  limit  is 
set  to  the  applicability  of  equation  97  by  the  corresponding  permeability 
limited  rate  of  transport  (see  6,4).  The  breakdown  of  equation  97  under  these 
conditions  is  obviously  due  to  the  neglect  of  viscous  forces  in  its  derivation. 

To  estimate  values  of  limiting  diffusion  rates  it  is  usually  sufficiently 
accurate  to  remember  that  T)^2  ®  P  being  total  pressure;  Dj^2  >  T 

being  absolute  temperature;  and  Dj^2  oc  s/  (1/mj^)  +  (l/m2) ,  nij^,  m2  being  the 
molecular  weights.  Thus,  assuming  Dj^2  about  0.1  cm.^/sec.,  t  =  0.1  cm., 

1  atmosphere  and  room  temperature,  n  =  2,  and  e/q  to  be  about  1/10,  then 

araps./sq.  cm.  Thus  for  equi-molar  counter  current  diffusion 
(f  “  1)  and  a  molar  fraction  of  reactant  near  1,  i  ~  1000  amps./sq,  cm. 

Ij 

This  figure  is,  of  course,  very  sensitive  to  t,  the  thickness  of  the  diffusion 
path  through  the  electrode.  For  air  at  the  same  conditions,  n  =  4,  Hj^q/N  -  1/5 
and  f  ■  1/5  therefore  i  2;  amps./sq.  cm.  An  increase  in  the  pressure  of  air 
does  not  produce  a  proportional  increase  in  the  limiting  current  because  although 
nj^Q  increases  proportionally  (and  f  remains  the  same),  ^-^2  inversel^y 
proportional  to  total  pressure  and  the  pressure  Increase  effe^  is  cancelled 
out . 

Knudsen  Diffusion 

When  the  mean  free  path  for  the  diffusing  molecule  is  greater  than  the 
pore  radius,  collisions  with  the  walls  provide  most  of  the  resistance  to 
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transport.  The  presence  of  other  molecules  will  thus  have  little  effect  on  the 
mass  transport  and  we  would  expect  the  diffusion  coefficient  to  be  independent 
of  total  pressure  and  partial  concentration,  but  dependent  on  the  size  of  the 
pores.  The  steady  state  of  motion  is  equation  91,  that  is, 

Sp/Sx  +  Sr/Sy  +  pv^v/Sx  =  0 

For  Knudsen  diffusion  where  the  viscous  forces  occur  at  the  surface,  v  is 
constant  with  y  except  at  the  surface  and  v  is  the  mean  mass  velocity  defined 
by  nv  =  J,  n  being  the  molar  concentration,  J  the  molar  flux.  Usually,  the 
velocities  of  flow  through  porous  media  are  low,  and  the  density  is  low,  so 
that  the  bulk  momentum  term  (or  kinetic  energy  term)  is  negligible  compared 
to  the  viscous  force. 

The  viscous  forces  are  due  to  the  change  of  momentum  of  a  molecule 
striking  the  wall.  The  loss  of  momentum,  in  the  direction  of  flow,  per 
sq.  cm.  per  second  gives  a  force 

F  =  (l/4)nCmvf 

where  f  is  the  accomodation  coefficient.  Substituting  J/n  for  v 

F  =■  (l/4)CmJf 

From  the  equation  of  motion,  for  a  tube 

trr^  dp/dx  *  -(1/4) JmCf 27rr 


J  *  -(2r/mCf) dp/dx 

(98) 

or 

J  =  D  dn/dx 

Iv 

(98a) 

where 

2;  (2/3)rC 

Thus  D  is  proportional  to  r,  to  ^/T  ,  and  to  1/  n/ m  but  independent  of 

Iv 

pressure  and  the  composition  of  other  gases  present.  Integrating  equation  98a, 
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11  =  (nFDj^/t)n^Q  (99) 

In  this  case^  an  increase  of  pr'essure  increases  n^^Q  but  does  not  decrease 
therefore  the  limiting  current  increases.  It  should  be  noted  that  Knudsen 
diffusion  and  Knudsen  permeability  are  the  same  physical  process;  when  only 
one  component  is  present. 

6. 3  Permeability 

When  a  single  gas  is  flowing,  in  laminar  flow,  through  a  tube 
of  diameter  d,  the  relation  between  flow  rate  and  pressure  drop  is  (18) 

=  (d^/32Ti)(p)(Ap/L)  +  (5Cd/3)(Ap/L)  (100) 

VIj  is  the  linear  velocity  at  Pj^,  T)  is  the  viscosity,  p  is  the  arithmetic  mean 
pressure  over  length  L,  Ap  is  the  pressure  drop  over  length  L,  8  is  a  constant 
close  to  one  and  C  is  the  mean  molecular  velocity.  The  first  term  in  the  RHS 
of  the  equation  derives  from  Poiseuille's  formula,  while  the  second  term  is 
the  slip  flow  term.  This  arises  from  Knudsen  flow  within  a  mean  free  path  of 
the  wall  which  changes  the  effective  viscosity.  As  d  becomes  small  the 
equation  reduces  to  the  Knudsen  equation,  with  an  appropriate  value  of  &  of 
8^.  As  d  becomes  large,  8  changes  somewhat  from  8^  but  the  slip  flow  term 
becomes  a  smaller  fraction  of  the  total  until  a  pure  Poiseuille's  equation 
is  obtained.  For  non-circular  flow  paths,  the  diameter  can  be  replaced  by 
4Mk,  where  M  is  the  mean  hydraulic  radius  and  k  is  a  constant  for  a  given  shape 
but  varies  with  different  geometries;  k  is  close  to  1.  Making  this  substitution 
and  rearranging, equation  100  goes  to 

p^Uj^  =  (M^/k^)(^/Tl)(Ap/L)  +  (M/k^(C)(Ap/L)  (100a) 

k^,  kj^  are  constants.  As  will  be  seen  later,  the  constants  characterizing  the 
geometry  and  pore  size  are  not  known  for  flow  through  porous  media  and  they 
have  to  be  determined  experimentally,  therefore  we  will  use  the  equation  in  the 
form 


Pj^Uj^  -  B^(p/T))(Ap/L)  +  KqC(Ap/L) 


(100b) 
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When  the  slip  flow  term  is  negligible,  is  the  permeability  coefficient 
(in  darcies  for  T]  in  centipoises  and  p  in  atmospheres)  and  replacing 
with  JRT,  p  with  NRT, 

J  =  B^(p/ti)(AN/L)  (100c) 

As  before,  J  is  the  net  molar  flux  density,  N  is  the  total  molar  concentration. 
For  a  one  component  system,  pure  reactant  transporting  to  the  reaction  site, 
for  example,  J  =  and  the  limiting  current  is  obtained  when  M  =  N.  Then 
p  *  P/2,  where  P  is  the  total  entrance  pressure  and,  using  t  in  place  of  L, 

i-  =  nF(B^/T])(P/2t)N  (101) 

Lt  O 

6.4.  Comparison  of  Diffusion  and  Permeability  Limitations. 

We  can  now  compare  the  predicted  limiting  currents  for  an 
almost  pure  reactant  transporting  through  an  electrode  with  no  gaseous 
product.  The  equations  are  97  and  101,  that  is, 

ij^  -  (nFDj2/t)N  ln(N/n2o) 

and 

=  (nFB^/tn)N(P/2) 

(n2Q  is  the  percentage  of  inert  in  the  feed  stream.)  The  smaller  of  these 
two  limiting  currents  will  be  the  best  approximation  to  the  actual  limiting 
current .  Comparing  the  two  values 

i  diffusion/i  permeability 
L  L 

=  110^^2  ln(N/n2o)/B^(P/2)  (102) 

2 

For  example,  let  us  consider  Dj^2  effective  to  be  about  0.01  cm.  /sec.,  T)  to 
be  about  100  micropoises,  P  Co  be  1  atmosphere  and  n2Q  to  be  1/10  per  cent. 

Then 


4 


1  dlffusion/i  permeability 
L  L 

=  (10‘^)(10"^)(2.3)(3)(2)/B^ 

'X.  10‘^/B  (102a) 

—  o 

Permeability  coefficients  of  0.1  to  0,40  have  been  reported  (22)  for  porous 
carbon  electrodes.  Thus  it  is  concluded  that  diffusion  will  almost  certainly 
be  the  mass  transport  limiting  mechanism  for  porous  electrodes. 

B^  is  related  to  the  viscous  area  and  flow  path  in  a  porous  solid.  In 
equation  100a,  is  the  true  velocity,  whereas  the  value  of  J  is  based  on 
the  external  face  area.  The  actual  open  area  for  transport  is  e  and  the 
actual  path  length  per  unit  thickness  is  q,  therefore  the  true  rate  of  flow 
in  the  solid  is  Jq/€.  Thus,  neglecting  slip  flow  and  assuming  ~  2,  equation 
100a  can  be  expressed  as 

J  =  (e/q^(M^/2)(p/Tl)M/L 

For  TJ  in  centipolses,  p  in  atmospheres  (1  atmosphere  -  10^  dynes/sq.  cm.) 

Bq  2;  (€/q^)(M^/2)(10®)  darcies  (103) 

M  is  the  effective  mean  hydraulic  radius  defined  by  equation  103.  It  is 
constant  for  a  given  material.  If  the  pore  size  of  the  material  does  not 
vary  widely,  M  can  be  approximated  by  an  overall  hydraulic  radius  defined  by 
total  cross-sectional  area  divided  by  wetted  perimeter  =  e/A,  A  being  as 

3 

before,  the  perimeter  area  per  unit  volume  of  electrode.  Thus  for  A  =  10 
2  3 

cm  /cm  ,  e  =  1/3,  q  =  2,  a,  1/2  darcy.  Again,  we  must  conclude  that 
equation  97  is  the  equation  which  will  be  applicable. 

In  general,  then,  we  may  assume  that  the  mass  flow  associated  with  net 
mass  transfer  can  be  accomplished  with  virtually  zero  total  pressure  gradient 
in  the  electrode  and  that  the  equation  to  be  solved  for  any  specific  set  of 
conditions  is  equation  95,  or  a  multicomponent  equivalent  of  this  equation 
(23). 
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6 . 5  Combined  Normal  and  Knudsen  Diffusion. 

Since  diffusion  is  the  predominant  factor  to  be  considered^ 
and  since  electrodes  containing  very  fine  pores  may  be  used  in  some  circum¬ 
stances^  it  is  worthwhile  to  investigate  the  combined  effects  of  the  two 
types  of  diffusion.  If  it  is  assumed  that  the  pores  of  different  sizes  are 
highly  interlinked  and  randomly  arranged,  then  it  can  be  assumed  that  the 
path  of  any  molecule  through  the  solid  will  involve  essentially  the  same 
resistances  to  flow  as  the  path  of  any  other  molecule.  Over  a  short  element 
dx  we  can  envisage  part  of  the  resistance  to  flow  being  caused  by  Knudsen 
effects,  giving  dn^^j^  concentration  drop  for  flow  J^. 

An  effective  Knudsen  diffusion  coefficient  can  be  defined  by 

Jl  -  .D^j^(dn^/dx)^ 

Similarly  an  effective  normal  diffusion  coefficient  can  be  defined  by 


+(n^/N)J 

The  total  pressure  drop  over  dx  is  dn^^j^  +  dn^/dx  •  dnj^jj)/dx. 

Therefore 


Thus, 


-dn^^/dx 


-  (n^/ND^2>'^ 


J.  =  -(l/(^  +  ^))dn  /dx  +  (n  /N)J(1/D  (^  +  ^)) 

^  *^12  '^Ik  ^  *^12  “^Ik 


(104) 


or. 


=  -Dj^  dn^/dx  +  (nj^/N)JB  (104a) 

For  »  0^2^  =  0^2  effective,  B  =  1.  For  D^2  '^Ik'  ^1  ~  ^Ik 

8  -»  0.  Thus  6  will  change  with  the  absolute  pressure  of  the  system.  If 
the  molar  flux  ration  J2^/J2  known,  then,  since  J  -  +  J2,  J  replaced 

in  terms  of  and  equation  104  can  be  Integrated  as  before  for  a  given  set  of 

circumstances . 
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(7)  FUTURE  WORK. 

This  report  lays  the  theoretical  grcound  work  for  a  number  of  inter¬ 
esting  investigations.  Examination  of  the  derivations  of  the  equations  shows 
the  assumptions  involved  and  it  is  then  relatively  easy  to  set  up  simple 
experimental  systems  which  approximate  as  neerly  as  possible  to  the  conditions 
required.  From  a  study  of  such  systems  it  shnould  be  possible  to  extend  the 
work  to  practical  fuel  cell  electrodes.  In  psract±cal  systems  it  is  inevitable 
that  more  complex  conditions  will  exist,  Forr  example,  the  flow  of  reactant  gas 
across  the  face  of  an  electrode  may  lead  to  c=oncerrtration  profiles  across  the 
face,  with  the  exit  gas  at  a  lower  concentration.  Even  if  the  flow  is  slow 
enough  that  free  space  diffusion  will  even  oust  th±s  concentration,  each  cell 
in  a  battery  of  cells  will  operate  with  a  different  feed  concentration.  Again, 
it  is  possible  that  high  output  cells  will  no»t operate  at  isothermal  conditions. 

Some  of  the  equations  derived  here  are  Welng  investigated  in  this  depart¬ 
ment,  but  not  enough  experimental  results  hav~e  been  obtained  to  date  to  warrant 
their  discussion.  Will  (12)  has  apparently  proved  the  validity  of  equation  63 
for  a  partially  immersed  platinum  wire,  with  hydrogen  fuel  and  sulphuric  acid 
electrolyte.  Allowing  for  what  is  certainly  a  misprint  in  the  abstract,  the 
form  of  equation  he  quotes  is  identical  to  th_at of  equation  63.  . 

It  is  to  be  expected  that  the  next  few  y  ears  will  see  the  development 
of  completely  adequate  theoretical  and  experirmental  descriptions  of  the  steady 
state  porous  electrode,  and  at  least  a  partia  1  description  of  unsteady  state 
situations . 
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Polarization,  volts  (b  >  0.026  volts) 
for  n  "  1,  room  temperature,  a  >  1/2. 


Polarization,  volts  (b  >  0.026  volts) 
for  n  •  1,  room  tenperature,  a  ■  1/2. 


Polarization,  volts  (b  >  0.026  volts) 
for  n  »  1,  room  temperature,  a  ■  1/2. 
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Figure  2D.  Polarization-Current  Relation  for  Equation  33. 
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Figure  5.  General  Shape  of  Equation  63. 
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Figure  9.  Decrease  of  Folarlzatlon  Versus  Tiaw  after  Current  Shut-off:  Irreversible.  (Equations  85^  86,  87,  88) 
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Figure  10.  Illustration  of  Concentration  Profiles 

In  Porous  Electrode  After  Current  Inter¬ 
ruption  (Reversible) . 


